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AN EVALUATION OF A CONSTRAINED TEST 
METHOD FOR OBTAINING FREE BODY RESPONSES 

By A. Berman, N. Giansante, F. D. Bartlett, Jr. 
Kaman Aerospace Corporation 


INTRODUCTION 


Dynamic testing, both full scale and model, is an essential 
step in predicting the response of aerospace vehicles to 
the conditions to which they will be subjected in flight. 
This testing is required for modal analysis, stability and 
control studies and loads analyses and is applied to design 
verification and modification studies. The actual in-flight 
boundary conditions, however, cannot be exactly duplicated 
on. the ground . 

In order to simulate the free-body boundary conditions of a 
vehicle in flight, the usual procedure has been to support 
the vehicle on a system which is relatively soft so that 
the "rigid body” frequencies (which should be zero) are low 
compared to the frequencies of the deformation modes of the 
structure. A commonly used technique (References 1, 2, and 
3) for launch vehicles consists of supporting the vehicle 
vertically on cables attached to its base. While tests 
conducted in such a manner seem to have given acceptable 
results, there are several disadvantages to this scheme. 

It is necessary to construct a tall structure capable of 
supporting the total weight of the vehicle. There is some 
uncertainty in the effects of the cable dynamics and non- 
linearities on the vehicle response (Reference 2). Various 
cable configurations have been known to give variation in 
test results (Reference 1). Certain new problems arise for 
vehicles which are not axisymmetric. When the center of 
gravity varies laterally under various fuel loads , the 
stabilization of such a vehicle on soft supports can become 
a major consideration. Vehicles which require testing in 
more than one attitude compound these difficulties. 

A procedure which could eliminate the effects of supports 
would be of significant benefit. It would not be necessary 
to use soft suspensions with the assumption that the 
interactions with the supporting structure are not signifi- 
cant. It would be possible to support the system being 
tested on a relatively stiff base, thus simplifying the 
problems of static stability and attitude variation. The 



design of supporting towers would be greatly simplified 
and the overall cost of testing would be reduced. It is 
essential, of course, that such a procedure be reliable, 
accurate, not overly sensitive to measurement errors, and 
applicable to real test conditions. 

There are several analytical methods which convert con- 
strained responses into free body responses . Typical 
methods are presented in References 4, 5, and 6. These 
methods, however, are suitable only for analytical pro- 
cedures where the response on infinitely rigid supports 
is known (or can be calculated) and where the mass matrix 
of the structure is available. Since such data is unmeasur- 
able in a test, these methods are not usable. 

The method which is examined here uses the measured forces 
of constraint to convert the measured structural responses 
to free body responses. The structure under test is con- 
sidered to be supported on real supports, but their specific 
characteristics are not required since only their measured 
reactions are used. The procedure uses only data which are 
actually measured, and no quantitative assumptions are used. 
The basis of this method was first discussed in Reference 7. 

The primary purpose of the work reported here is to evaluate 
the suitability of the method for use under realistic con- 
ditions and for representative aerospace vehicles. Addi- 
tional purposes are to establish guidelines for usage of 
the method: and to provide computer software capable of 
analyzing actual test data. 
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DESCRIPTION OF THE THEORY 

> 

BASIC CONCEPT 

Consider a constrained structure which is being shaken by 
a known force and assume that the reaction forces at the 
supports are known. The structure responds precisely as 
if it were a free body being. simultaneously subjected to 
the actual applied forces and to the forces of constraint. 
Thus,. a shake test in which the constraining forces are , . 
measured gives direct information about the free body re- 
sponse of the structure when acted upon by several forces.. 

As will be seen below,. it is possible to convert . information 
of this type into the free body response of the system when 
subjected to only one force at a time. This is what; is 
needed to determine resonance data and to predict the 
effects of arbitrary loads. 

ANALYSIS 

The analysis assumes that the structure being studied is 
linear at each frequency. That is, the response is pro- 
portional to the force and the principle of superposition 
holds. For sinusoidal forces at a frequency, w, applied 
at M points on a structure, there is defined a vector f^ , 
which represents the complex amplitude of applied force 
at each of the points. Similarly, y^ , is a vector rep- 
resenting the complex amplitudes of the deflection at each 
of N points resulting from the force, f^. There .is no 
necessity for the force points (represented in -T^^to 
coincide with the response points (represented in y^) . Y 
is the complex NxM displacement mobility matrix representing 
the relationship between the forces and responses. F^, y^, 

Y and the relationship between them are written: 

3y x 3y x 3y x 

3y 2 : 

3 f i 

^N 9 *N 

9f l 9 *M 
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and 


y k = Yf k (1 » 

The displacement is used only for illustration, exactly the 
same relationships hold for velocity and acceleration. The 
displacements can also, with no change in the analysis, rep- 
resent displacements or rotations in two or three directions 
at one geometrical point by allowing one element in each 
vector for each of these generalized displacements. Similar 
considerations apply to the forces (or moments) . Note that 
there is no necessity for Y to be square, it will contain 
one row for each displacement measured and one column for 
each point at which a force is applied and, as will be seen 
below, one column for each constraint. 

Consider, now, a matrix, F, containing several applied load 
vectors and a matrix, y, containing the corresponding de- 
flections , as follows : 


F ■ t? i 7 2 ** •] y = [Yi Y 2 * • •] 

and then 


y = YF (2) 

If F is a nonsingular matrix, then the desired result, the 
response of particular points to single forces, may be 
written 

Y = yF -1 (3) 

where both y and F are measured. When the "actual" applied 
loads only are included in F, then Y is the mobility of the 
structure as tested - i.e. on the actual supports. If F 
includes any of the forces of constraint, then Y is the 
mobility of the structure with those constraints removed. 

If F includes all the forces of constraint, then Y is the 
mobility of the free body. ~ 

As stated above, F must be nonsingular and thus have an 
inverse. If there are M forces to be considered (including 
the forces at the constraints) then M sets of forces, 
must be applied and all of these vectors must be in- 
dependent. There are at least two ways that this may be 
done: (1) by applying an external force at each constraint 

or (2) by varying the constraints. 
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FORCES AT CONSTRAINTS 


If an exciting force is applied at the k-th constraint, the 
force vector will be of the form 



(4) 


where the r's are measured forces of constraint and f is the 
applied force. The force vectors obtained by applying forces 
to the structure at each of the constraint points will 
ordinarily be independent of each other and the force matrix 
will, therefore, be nonsingular. In the event the force 
matrix is ill-conditioned, the location or characteristics of 
the constraints may be altered to yield a well behaved force 
matrix. The matrix of forces at the supports including the 
applied forces will be called F s . 


At the same time that these forces are measured, the dis- 
placements are measured at the points of interest on the 
structure and one column of y is formed for each column of 

F . Then , as above 
s 


y = yF s _1 (5) 

where Y represents the deflection of each point of interest 
due to each of the loads applied to the structure (at the 
supports). This is the free body mobility matrix. This 
procedure must be carried out over the frequency range of 
interest. 

If it is desired to find the response due to forces applied 
to the structure at points other than the support points, 
then the structure must be shaken at these points in addition 
and the forces at the constraints must be recorded. If Fc^ 
is a matrix representing the forces of constraint for each 
non-constraint point of excitation, then the F matrix becomes 
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( 6 ) 


where unit forces are applied. The inverse of this matrix 
involves little more than inverting F and is given by 

d 



(7) 


In practice it will often be just as convenient to invert 
the matrix of Equation (6) directly. Equation (7) illus- 
trates that no numerical complications are introduced by 
forces at additional points. 

VARIED CONSTRAINTS 


Any means of varying the constraint forces such that the F 
matrix is nonsingular will work. Applying a force at each 
of the constraints was just discussed. Another possible 
method is to vary the constraints themselves such that the 
force vectors are independent. 

If the structure is supported redundantly, then a procedure 
which would work is to shake at only one constraint and 
remove one constraint at a time resulting in an F matrix of 
the following form (where the force is applied at station 1) 


f+r^ f+r^ f+r^ f+0 f+r^... 


F = 


( 8 ) 


where the first column represents the measured loads when 
all the supports are used, the second represents the loads 
with constraint number 3 removed, etc. 
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The same effect can be achieved by varying some parameters 
e.g. the stiffness, of each constraint one at a time. Thi 
would eliminate the need for redundant supports and reduce 
the amount of data required. 

The other considerations are similar to the previous 
method . 



APPLICATION CONSIDERATIONS 


The method has attributes which make it an especially 
attractive candidate for practical application. These 
include the use of only measured data and the lack of 
quantitative assumptions . There are, however, as in all 
procedures , certain considerations involved in planning 
an efficient and accurate application of the method. 

NUMBER OF CONSTRAINTS 

At each frequency it is necessary to conduct one test for 
each constraint, thus it is desirable to keep this number 
to a minimum. While it is possible to constrain all rigid 
body motions with six constraining forces, there is no 
necessity for such complete constraint, however. During 
the design of a test, consideration should be given to test 
configurations which allow freedom of motion, e.g. in the 
horizontal plane and around the vertical axis. In this 
case, it would be necessary to shake vertically at each 
support and measure each of the vertical forces of con- 
straint. In addition, any other shaker position or 
orientation could be used while the three vertical forces 
were measured. 

For the design of a specific test, it is necessary to 
evaluate the cost of eliminating constraints compared to 
the reduced testing required. 

SUPPORT CHARACTERISTICS : 

In a theoretical sense, the characteristics of the supports 
are immaterial. These characteristics, however, do affect 
the magnitudes of the forces and displacements which will 
be measured. The accuracy of the various transducers 
(accelerometers, load cells, etc) depends on the magnitude 
and frequency of the quantity being measured. Thus, for 
the most reliable results , the supports should be designed 
and the transducers selected so as to operate in their 
most accurate region. It is not expected that this will 
be an extremely stringent requirement, but care must be 
exercised as in planning any test, to insure that the 
results be meaningful. 
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In the analysis, the constraint forces are considered to 
be sinusoidal and at the same frequency as the excitation 
force. For real supports, it is not uncommon for components 
of other frequencies to pollute the response. It is 
implicitly assumed here that any such components have been 
removed from the data by electronic or digital filtering 
or by Fourier analysis. Obviously, highly nonlinear supports 
will increase the uncertainty in the data and should be 
avoided. . 

ERRORS 


Any correctly derived analytical procedure will work when 
the input data is exact. When such methods are applied to 
measured data the effect of the unavoidable experimental • 
errors is critical in the evaluation of the practicality 
of a particular technique. This process uses a .measured 
force matrix containing errors, inverts this matrix and 
multiplies by another measured matrix also containing 
errors. The behavior of these errors will determine 
whether the method, is economically feasible. The 
expected accuracy of the final results . compared to 
the accuracy of an alternate method is an important 
consideration. ; 

TEST SIMULATION 

A computer simulation of alternative test configurations 
can be an extremely useful tool in the preliminary design 
of any test. Such a simulation applied to this method 
using an approximate analytical model of the vehicle can 
be used to determine the sensitivity to error and the ex- 
pected accuracy of .the results of (the various arrangements . 
considered. It should include realistic experimental 
errors, approximate constraint characteristics, and vary 
the frequency over the range of interest. - 



SIMULATION SOFTWARE 


The computer software developed under this study has been 
designed to serve several purposes . The programs can be 
used to perform a general evaluation of the feasibility of 
these techniques. This is the particular application made 
in the work reported here. 

In addition, the program has been designed so as to be 
capable of analyzing actual test data. This facility was 
accomplished by dividing the software into two separate 
logical entities. The first develops simulated test data, 
the second analyzes the test data which can be either 
simulated of actual. 

A third application of the software is for the planning of 
an actual test. This usage requires approximate analytical 
models of the structure to be tested and the constraints. 
The programs can be used to optimize the locations and 
general characteristics of the constraints and to estimate 
the accuracy of the resulting data. 

COMPUTER PROGRAM ORGANIZATION 

The portion of the software used to develop the simulated 
test data is itself divided into two separate programs for 
efficiency reasons. The first program takes as input the 
mass and stiffness matrices of the structure being modeled. 
The structural damping coefficient is also specified. This 
quantity is allowed to vary over the structure. If g^ is 
the damping coefficient at coordinate i then the diagonal 
elements of the imaginary damping matrices are gi ven by 
9i K ii an< ^ the off-diagonal elements are given by /g^g • no- 
where and K^j are diagonal and off-diagonal 
elements of the stiffness matrix. In addition, and as 
separate inputs for convenience, the support characteris- 
tics are entered. These characteristics include the 
stiffness, structural damping, and/or viscous damping 
coefficient. 

These quantities, of course, would be unknown during an 
actual test but they must be estimated for purposes of 
simulation. The program computes exact values of the 
constrained system mobility matrices at a set of specified 
frequencies and writes this data on tape. 
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The second program uses the exact data computed by the first 
program and introduces specified errors (see below) to yield 
simulated acceleration and force measurements. The simula- 
ted test data is written on a tape in a format which is 
compatible with actual test data. The exact data developed 
in the first program may be reused. with different errors 
or forces. Because the first program is by far the most 
time consuming, this feature improves the efficiency of 
the simulation process. 

The third program uses as input, either simulated data (as 
developed above) or actual (but pre-processed) test data. 
This program reduces the measured constrained accelerations 
and forces of constraint to yield the free body responses 
of the system. 

Detailed listings of the programs, description of input 
formats and the various options available are given in the 
appendix. 

SIMULATED ERRORS 

In order to simulate test data for the evaluation of a 
numerical procedure, it is essential that the simulated 
data contain realistic errors. 

The program has provisions for several types of errors to 
be simultaneously applied to each of the simulated measured 
accelerations and the simulated measured forces. These 
errors include: (1) a uniformly distributed random per- 

centage error on amplitude between specified limits; (2) a 
uniformly distributed random phase angle error between 
specified limits; (3) a constant specified percentage bias 
error; (4) a uniformly distributed random amplitude error 
between limits (simulating system noise) . 


The simulated measurement errors which are estimated to be 
representative of the current state-of-the-art and applied 
in the computations were generally as follows: +5% random 

error on amplitude of acceleration and force; +5^ random 
phase error on acceleration and force; a random +2.54 cm/sec 
(+1 in/sec 2 ) and a random .01 radian/sec 2 on translational and 
rotational acceleration measurements respectively. However, 
several computer experiments were conducted using error 
values ranging to +15% random and +15% random phase error 
to test the sensitivity of the method to error. 
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MODELS ANALYZED 


The analytical models analyzed were representations of a 
1/15 scale dynamic model of a space shuttle configuration 
(described in Reference 8) . The various stiffness and 
consistent mass matrix formulations corresponding to each 
of the structures considered were obtained using NASTRAN 
and were supplied by NASA. In addition, structural damping 
coefficients of 3% and .5% were used for the orbiter and 
booster fuselages respectively. The basic models studied 
are illustrated in Figures 1 and 2. 

Several models were investigated including: an eight 

coordinate orbiter fuselage limited to transverse motion; 
a two-dimensional or 16-degree-of-f reedom representation 
of the orbiter fuselage with each of the eight coordinates 
possessing a transverse and rotational degree of freedom; 
a three-dimensional model of the orbiter fuselage; with each 
point having six degrees of freedom, yielding a total of 
48 degrees of freedom; a coupled orbiter and booster model 
consisting of 18 coordinates with each point having trans- 
verse and rotational degrees of freedom; a one-dimensional 
model of the booster fuselage having 10 degrees of freedom. 



DISCUSSION; OF RESULTS 


Computer simulations of dynamic tests were conducted to 
obtain the .free body response of the aforementioned models. 

A variation" of parameters was considered including magnitude 
and location of constraints, structural and viscous damping 
of the constraints and experimental error. The simulated 
tests were carried out at various frequencies over a spec- 
trum of interest peculiar to each model studied. 


Except when otherwise noted, the errors used were as follows 
uniform random amplitude error of +5%, uniform random phase 
error of +5°. uniform random absolute amplitude error of 
+2.54 cm/sec 2 (+1 in/s ec 2 ) on translational accelerations 
and .01 rad/sec 2 ^ on rotational accelerations. Also, the 
corresponding force random errors used were +5% and +5°. 

Each simulated test was run with force levels of 4.45, 22.24 
44.48, 111.21N (1, 5, 10, 25 pounds). Each of these was 
repeated with different random sequences. Thus, every range 
of data shown in the following plots represents the extreme 
values obtained in 8 simulated tests. The predominance of 
the data was obtained using the procedure of shaking at 
the constraints. 


For convenience the plots of acceleration response per unit 
force are shown without dimensions since only relative 
amplitudes are of interest here. However, for reference, 
unity represents .0057 m/sec 2 /N (1 in/sec 2 /lb) 


A large number of simulations were computed. Only a small 
portion of the data for each of approximately 20% of the 
cases run is presented in this report. The conditions 
shown were selected to illustrate typical results. The 
remaining data, which is available (and easily duplicated) , 
would be purely repetitive and add no new information. 


TWO-DIMENSIONAL ORBITER 


The two-dimensional orbiter is a relatively simple model 
and much of the early exploratory work was done with it. 
Some typical data is presented here. Figure 3 presents 
for reference the exact normalized acceleration amplitude 
frequency response and the associated phase angle of the 
two-dimensional orbiter in free and constrained configura- 
tions. For the constrained system, the supports consisted 
of both 1.75 x 10°N/m (10,000 lb/in.) translational springs 
and 1129 m-N/rad (10,000 in-lb/rad) rotational springs ac 
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stations .69m and 1.66m, with zero structural and viscous 
damping. Application of the theory to the simulated test 
constrained response yielded the free body characteristics 
of the structure. The calculated natural frequencies for 
the range of frequencies shown are in excellent agreement 
with those supplied by NASA. 

Figure 4 illustrates the effect of including measurement 
error in the simulated measured accelerations and the simu- 
lated measured forces. The ranges of the free body response 
with errors included are shown superimposed on the exact or 
zero error free body response. The error bounds indicated 
at the various frequencies were the maximum and minimum 
values obtained from several computer runs with the same 
nominal error, but with different random number seeds. The 
consistency of results from individual simulated tests 
indicates the method is relatively insensitive to the level 
of nominal errors applied. The natural frequencies are 
accurately defined and the general shape of the response 
curve is retained even at bounds of the errors. In Figure 5 
a typical simulated test frequency response is shown com- 
pared to the exact response. It is apparent that the exact 
response and the response deduced from test are in excellent 
agreement . 

The amplitude of the responses of the structure at 92 Hz, 
which is approximately at the first natural frequency of 
the free system, due to a force at station 1.66m is 
presented in Figure 6 . The constrained response and the 
free body response are shown for the exact conditions. 

The free body responses including the effect of measurement 
error are indicated in the figure as ranges of values . 

THREE-DIMENSIONAL ORB I TER 

The three-dimensional orbiter model represents a real struc- 
ture with six degrees of freedom. The data presented repre- 
sents the structure on six constraints and illustrates the 
effect of the stiffness of the supports. Figure 7 presents 
the constrained normalized acceleration amplitude frequency 
response and the accompanying phase angle for the three- 
dimensional orbiter model restrained with 8.75 x 10 5 N/m 
(5000 lb/in) springs in the transverse, lateral and 
longitudinal directions at station .69m and 1.66m. The 
exact free body responses are also shown. The result of 
applying, simultaneously, the several types of errors to 
both the simulated measured accelerations and the simulated 
measured force is indicated in Figure 8. Error bounds at a 
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particular frequency were obtained by several computer 
experiments at the same nominal error levels but varying 
the random number seed and applied force. For the error 
values considered the free body response does not vary 
significantly from the exact values. Figures 9 and 10 
present similar data except that the free body response 
was obtained from the constrained configuration with 1.75 x 
10 6 N/m (10,000 Ib/in. ) and 8.76 x 10 6 N/m (50,000 lb/in.) 
spring rates , respectively. A comparison of Figures 8, 9 
and 10 reveals that, as would be expected, the free body 
response deduced from the constrained structure, for the 
zero error condition, are identical regardless of constraint 
configuration employed. The figures also illustrate that 
the larger ranges in results occur in the response regions 
of extreme slope. However, at these points the error bounds 
are incidental in defining the actual response curve. 


The data shown in Figure 11 is the transfer response, the 
transverse acceleration at station .31m due to force ex- 
citation at station .69m. The exact response of the struc- 
ture at 90 Hz, which is slightly below the first natural 
frequency of the free system, due to a force at station 1.66m 
is presented in Figure 12 for the free condition and for 
the restrained structure with constraints at station ,69m 
and 1.66m. The results of the calculated free body responses 
are shown as ranges of values. Figure 13 illustrates a 
typical free body response converted from constrained data 
including measurement error compared to the exact free body 
response. Figures 14 and 15 present constrained and free 
body responses similar to Figure 13, however, the constraint 
spring rates were 1.75 x lO^N/m (10,000 lb/in.) and 8.75 x 
lO^N/m (50,000 lb/in.) respectively. A comparison of 
Figures 12, 14 and 15 illustrate the significant variation 
in free body responses obtained from the simulated measured 
test data, with error, for the different support conditions 
considered. The results of using the lower spring rate are 
more consistent and vary least from the zero error response 
although the same, nominal errors were used in the simulated 
measured data obtained from each support system. In regard 
to Figure 15, it should be noted that in the frequency 
response plot, Figure 10,' the frequency always resulted 
in values below the exact. This is the only frequency at 
which this condition occurred. 
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The results of the deduced free body responses are 
presented as ranges of values for ' each ' particular frequency 
at which the simulated tests were conducted. There is no 
apparent deterioration in the results due to this constraint 
condition. 

To further examine the effect of experimental error on the 
analysis, simulated tests were performed using error levels 
ranging to +15% random error on amplitude of accelerations 
and forces and +15° random phase error on the same param- 
eters. Only one simulated test was run for each level of 
error at the indicated frequencies and the results are shown 
on Figure 25. Also presented in the figure are the ranges 
of values for the simulated test with nominal measurement 
error levels as described previously. 

BOOSTER LONGITUDINAL RESPONSE 

This model is an extremely simple one but does correspond 
to actual tests. Here the data is used to illustrate that 
damping in the supports has no ill effect on the results. 
Figure 26 presents the frequency response of the ten-degree- 
of-freedom booster fuselage in the free condition. The 
exact amplitude and the exact phase angle response are shown 
with the results including error in the simulated measured 
data given as ranges of values. The free body response was 
extracted from simulated measured constrained data which 
incorporated zero structural damping in the restraint system. 
It is apparent that the free body response is insensitive to 
the level of applied error, at least for the conditions 
investigated. Figure 27 presents data similar to that shown 
in Figure 26, however, the free body response was derived 
from simulated measured constraint data with 5% structural 
damping in the supports. Further computer experiments were 
conducted using 20% structural damping in the restraints. 

The resulting error bounds superimposed on the exact free 
body ‘response are effectively the same as those shown in 
Figure 27, therefore, the range of values for the 20% 
structural damping condition are not presented. Based' on 
the number of conditions analyzed and the levels of con- 
straint structural damping used, it appears the free body 
response of the structure is insensitive to the level of 
structural damping in the constraint system. 

The exact free body response and the constrained response 
from which it was deduced are given in Figure 28. Also 
shown are the error bounds resulting from the computer 
experiments in which simulated measurement errors were 
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considered. The responses, due to force excitation at 
station 3.41m, were calculated at 445 Hz which is the 
approximate second elastic natural frequency. 

IGNORED CONSTRAINTS 

The forces at each constraint must be measured in order to 
obtain the response of the completely free body. It is also 
necessary to force at each of the constraints when this 
procedure is being used. In practice, there may actually be 
constraints in directions where none had been planned. For 
example, a pin joint under load may exhibit a torsional 
restraint. The question arises as to whether this method 
will deteriorate if some of these forces are ignored, either 
accidentally or intentionally. 

It is necessary to reexamine the theory to answer this 
question. It has been shown that the responses deduced 
by this method are those for the structure with those 
constraints removed whose reaction forces are included in 
the force matrix. Thus, if constraint forces exist and 
are ignored, the computed responses are those of the 
structure on those constraints only. When the ignored 
forces are small, they will have a small influence on the 
free response. This situation is quite analogous to the 
soft suspension concept for free body testing (where, in 
actuality, the forces of constraint are ignored) except 
that here , the remaining constraints are not required to 
provide static stability to the structure being tested. 

Thus, when the forces of a constraint are small enough to 
have a negligible effect on the free body response they 
may be ignored. 

Several simulated tests of this condition were run using 
the two-dimensional orbiter model constrained torsionally 
as well as translationally (see Figure 3). Program 3 has 
a capability for masking selected channels of data on the 
input tape and this facility was used to ignore the tor- 
sional forces produced at the constraints. In these cases, 
the differences between the computed responses and the 
exact free body responses were too small to illustrate. 

These tests confirm the theoretical conclusion stated in 
the previous paragraph. 
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NUMBERS OF CONSTRAINTS 


Some of the simulated tests illustrated had more constraints 
than necessary. Some allowed freedom in one or more 
directions . 

The two-dimensional orbiter data presented here (Figures 
3-6) actually had four constraints where only two were 
needed for full constraint of the system. Other tests not 
illustrated used only the two linear springs. There 
appeared to be no loss in validity through the use of 
redundant constraints. Redundant constraints, of course, 
increase the amount of testing required. 

The number of constraints also appear to have no particular 
degrading effect on the data. The three-dimensional orbiter 
data (Figures 7-16) used six constraints with results as 
good as data obtained from lesser numbers of supports. 

Allowing unconstrained motion in one or more directions 
also was found to have no noticeable effect on the data. 

The three-dimensional orbiter tests did not constrain the 
structure in roll and no loss in accuracy was observed. 

From a theoretical point of view none of these conditions 
was expected to affect the final accuracy and the data 
obtained verified this conclusion. 

VARIATION OF CONSTRAINTS 

Figure 29 illustrates the result of applying the variation 
of constraint technique described in the theoretical de- 
velopment to the two-dimensional orbiter fuselage model. 

The data was achieved by varying the constraint spring rate 
while maintaining the same force excitation on the struc- 
ture. Initially, the structure was restrained with 1.75 x 
10°N/m (10,000 lb/in.) springs at station .69m and 1.66m 
with force excitation applied at station 1.66m. Subsequently 
the spring rate at station .69m was changed to 3.5 x lO^N/m 
(20,000 lb/in.) while the spring rate and force level at 
station 1.66m were held constant. Figure 29 presents the 
free body exact acceleration response of station 1.66m to 
force excitation at the same station. Also shown are the 
deduced ranges of values obtained from the simulated con- 
strained tests applying error to the measured forces and 
accelerations . 
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The results shown on Figure 29 compare favorably with those 
presented in Figure 4 which were obtained by the method of 
force application at each of the constraints while the 
support characteristics remained constant. On the basis 
of sensitivity to measurement error there appears to be no 
difference in the results obtained from the two methods. 


i 
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CONCLUSIONS AND RECOMMENDATIONS 


GENERAL 


The method of obtaining free-body responses through con- 
strained testing was investigated for practical applicability. 
A large number of simulations of shake tests were carried 
out using several representative analytical models of space 
shuttle components and systems. Constraint locations and 
characteristics were varied. The procedure appeared to work 
quite well (subject to considerations discussed below) and 
its application to actual testing is recommended. 

ACCURACY 

The accuracy of the deduced free-body responses is believed 
to be as good as or better than other procedures using soft 
suspension systems where the constraining forces are 
assumed small and are ignored. The random errors assumed 
in the test, simulation are believed to be within the present 
state-of-the-art of acceleration and force measuring systems. 

For all the conditions presented, the free-body amplitude 
and phase data obtained versus frequency adequately defined 
the response curves including resonances and antiresonances. 

SENSITIVITY 


The data obtained suggest three general conditions under 
which the computed free-body response may be especially 
sensitive to measurement errors: 

(1) , At a frequency near a free body resonance which 
is not near a constrained resonance, larger errors occur. 
This is not unexpected since small responses must be con- 
verted into large ones and errors can be expected to be 
amplified. In general, the data has indicated that even 
with relatively large uncertainties, the frequency response 
is still well defined because it is nearly vertical in 
this region. This situation will only occur at a small 
number of lower modes, because at higher frequencies the 
constrained and free responses become closer to each other. 

(2) Constraints which are in close proximity tend to 
increase the error sensitivity. This effect is due to the 
fact that constraints which act in the same direction and 
are very close to each other will tend to have nearly the 
same forces of constraint. This causes the corresponding 
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columns of the force matrix to be nearly equal and the 
matrix to be ill-conditioned. In general it is good 
policy to keep the constraints well separated when 
possible. However, this is not always possible, for 
example, when a space vehicle is in launch attitude and 
it is resting:^ on several supports at its base. Under 
these conditions the supports should not be considered 
to be independent but should be treated as a single support 
with up to three perpendicular force constraints and up 
to three moment constraints. This treatment will eliminate 
the ill-conditioning of the force matrix. There is one 
condition which should be avoided. This is when the 
supports are not close enough to be treated as a single 
support, as above, and yet not far enough separated to 
yield a well-conditioned force matrix. This condition is. 
expected to be rare. 

(3) The third situation which appears to aggravate 
the error sensitivity- is when the constraints are excessively 
stiff. This, also, is not unexpected. There are two reasons 
for this effect. First, the stiffer the supports, the 
further the actual response will be from the free-body 
response and the more the data has to be modified. Second, 
the stiffer the supports, the smaller the response of the 
structure for the same force and the greater the error to 
be expected in the raw data. This effect will show up in 
the data obtained in these simulations because of the noise 
type error' included in addition to the percentage amplitude 
and the phase errors used. 

The use of very stiff supports is to be discouraged, in 
general, if for no other reason than to reduce the exciter 
force requirements. In general, the supports should not 
be so soft that there is only a small difference between 
free and constrained responses because of the problems 
associated with static stability, etc, as discussed in the 
Introduction. On the other hand, they should not be so 
stiff that excitation of the structure will be a problem. 
There ’.is a very large middle ground where the techniques 
studied here are especially applicable. 

PRE-TEST SIMULATION 

Prior to applying these procedures in an actual test, it is 
recommended that an approximate simulation be carried out. 

The computer programs supplied are perfectly suited for 
this purpose. The input required is • a simple model (a 
highly reduced NASTRAN model, for example) and very 
approximate characteristics of the supporting system under ' 
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consideration. The simulation will determine whether any 
of the conditions mentioned above will be approached and 
will allow the consideration of alternate supports or 
support locations which will improve the validity of the 
data. 

Simulation is not essential and confidence may be had in 
the results if the test is set up with the sensitivity 
considerations, above, in mind. However, such a simulation 
would be an inexpensive precaution and would probably 
result in more valid data. 

POST-TEST SENSITIVITY ANALYSIS 

The program supplied will convert test data into free-body 
data after the test data has been calibrated and filtered 
and properly formulated on a tape. The particular tape 
format used was selected to make this preprocessing as 
simple as possible. 

The addition of a feature to this program is recommended 
which will establish confidence limits on the free-body 
data. The procedure would be as follows: (1) the user 

supplies his best estimates as to the ranges of errors 
contained in the input data; (2) the program then auto- 
matically and randomly varies each element of input within 
the specified limits and produces a distribution of 
variances in the deduced free-body responses in a form 
similar to that of Figure 16. This provides an estimate 
of the possible error in the free-body response. 

IMPLEMENTATION 


The primary method considered in this report calls for 
shaking the system at each constraint. The literal 
implementation of this system may entail some problems 
of convenience. Unless multiple shakers are to be used, 
probably the most feasible method would be to install a 
shaker at one support, perform a frequency sweep, then ; 
move it to another support and perform another frequency 
sweep until all the required data was obtained. There are 
two disadvantages to this procedure: (1) the free-body 

responses cannot be obtained until all the testing has been 
completed; and , (2) care must be exercised to insure that 
the data is taken at precisely the same frequencies. 
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In order to implement this procedure , the supports must' be 
designed to allow for application of the shaker. It is 
suggested that only linear {and no torsional) constraints 
be used since moment applications would not be as convenient 
to apply or measure. A type of constraint worth considering 
is the cantilever flexible beam as described in Reference 9. 

It would be possible to mount a shaker directly beneath the 
support. Because the beam would not have the extreme 
flexibility requirements as in the reference , it does not 
have to be cantilevered but may be supported at two or more 
points. Another advantage to the beam support is that it 
is possible to easily move the shaker from one position to 
another by sliding it along the beam. 

If several shakers are available it is possible to obtain 
all the data needed at one frequency before proceeding to 
the next frequency. This allows the possibility of a simple 
real time data reduction system yielding on site visualiza- 
tion of the free-body responses. This possibility includes 
the conveniences of "free-body" testing with the stability 
and adaptability of constrained testing. 

It is recommended that a simple implementation of this 
technique be tested in practice and that an evaluation of 
more sophisticated implementation procedures be carried 
out for future use. 

i 
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Figure 2. Orbi ter-Boos ter Combination Model 
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Figure 3. Frequency Response 
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Figure 4. Effect of Error on Frequency Response 
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Figure 7. Frequency Response 
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Figure 8. Effect of Error on Frequency Response, 
8.75 x 10^N/m Restraints 
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Figure 9. Effect of Error on Frequency Response, 
1.75 x 10 6 N/m Restraints 
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Figure 11. Transfer Frequency Response 
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Figure 12. Transverse Bending Response 
at 90 Hz, 8.75 x io 5 N/m 
Restraints 
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Figure 14. Transverse Bending Response 
at 90 Hz, 1.75 x 10 b N/m 
Restraints 



40 









TRANSVERSE ACCELERATION 


200 


100 


0 


-100 


STA (p.) 



THREE-DIMENSIONAL ORBITER 


Figure 15. 


Transverse Bending Response 
at 90 Hz, 8.75 x 10 6 N/m 
Restraints 


80 Hz hi ■ 90 Hi ' u e 140 Hz 



g 


E-* « 
Z W 

Bo 

cd w 

W H 

ws. 




o o o o o o 

o as so ««r in 

rH 

H0HH3 Q3XVDIQNI N1HXIM SXN3W313 JO 1N3DH3J 


42 


Figure 16. Error Distribution for Free Response 
3D Orbiter Model , Constrained at 
Stations .69 and 1.66 Meters '• 
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Figure 20. Transfer Frequency Response 
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Figure 26. Effect of Error on Frequency Response 


52 










Q 

W 

Z 

M 


£ 



NOiiLYHanaDDY avNiannoNoa 


54 


ONE-DIMENSIONAL BOOSTER 
Figure 28. Longitudinal Response at 445 Hz 








COMPUTER PROGRAM DESCRIPTION 


The digital computer program was designed to test and evalu- 
ate the concepts presented in the theoretical development. 
The program generates simulated test data which is as re- 
alistic as possible and operates on this data to yield the 
free body response of the structure, The program was 
written using CDC FORTRAN 2.3 language and can be run on 
the CDC 6400 and 6600 computer. 

The program is limited presently to sixty degrees of 
freedom. However, with slight modification this restriction 
can be removed depending upon the storage capacity of the 
computer used for program implementation. 
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GENERAL NOTES 


All integer variables must be right justified with no 
decimal point. 

• * 

Tape Assignments 

IT1 = 1 contains complex mobility matrices at specified 
frequencies . Tape IT1 is used as input to 
PROGRAM K2LRC. ‘ 

IT2 = 2 contains constrained complex acceleration matrices 
and simulated measured force data. Tape IT2 is 

used as input to PROGRAM K3LRC. 

» 

All input data must be in consistent units, frequency in 
Hertz. 

Card reader used as input devices . 

Printer used as output device. 
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PROGRAM K1LRC 
INPUT 


Card 


1 Columns 1 


2 

3 

4 

5 

6 


7 

8 


Parameter 


Definition 


IC1 = 0 Normal Program Execution 

Read Mass, Stiffness and 
Constraint Data 
= 1. Read New Constraint Data 
After Initial Run 
= 9 End Program 


IC2 

IC3 


1 Not Used 


IC4 


IC5 


1 Diagonal Mass Matrix 
= 1 Full Mass Matrix 


IC6 


IC7 

IC8 


= 0 Constant Structural 
Damping 

5 * 0 Variable Structural 
Damping 

| Not Used 


HEADC Program Heading 


Card 2 Columns 1-10 ND Number of Degrees of Freedom 


Card(s) 3 


11-20 GG Constant Structural Damping 

AM Mass Matrix 

If IC5 ? 1 Input Diagonal 
Mass Matrix, Cards Contain 
Diagonal Elements in 
Sequence, Format (8E10.0) 

If IC5 = 1 Full Symmetric 
Mass Matrix Input Lower 
Triangular Portion Starting 
Each Row on a New Card. 

Use as Many Cards as Nec- 
essary, Ending With Diagonal 
Element Format (8E10.0) 
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Parameter 


Definition 


Card(s) 4 AK 


Card(s) 5 G 


Card 

6 

Columns 

1-10 

NCON 

Card (s) 

7 

Columns 

1-10 

NCOR 




11-20 

AKR 




21-30 

GR 




31-40 

DAMP 

Repeat 

Card 7, One 

i Card 

for Each 

Card 

8 

Columns 

1-10 

NF 

Card (s) 

9 

- 


HZ 


Stiffness Matrix. Input 
Lower. Triangular Matrix 
Similar to Procedure for 
Mass Matrix Format (8E10.0) 

If IC6 ^ 0 Input Variable 
Structural Damping Vector, 
One Element for Each Degree 
of Freedom Format (8E10.0) 

Number of Constraints 

Coordinate Number at 
Constraint 

Constraint Spring Rate 

Constraint Structural 
Damping 

Constraint Viscous Damping 
Constraint 

Number of Frequencies at 
Which Mobility Matrices 
Will be Calculated 

Frequency Values at Which 
Mobility Matrices Will be 
Calculated Format (8E10.0) 
Ten Columns Per Value, 

Eight Values Per Card 
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PROGRAM KlLRC - SUBROUTINES 


SYM Forms symmetric matrix from lower triangular matrix. 
Uses object time dimensions. 

Used to form symmetric mass and stiffness matrices 
from corresponding lower triangular matrices. NRA 
is the dimensioned number of rows of the matrix and 
N is the order of the square matrix being considered. 

M0UT2 Used to print real matrices. 

Uses object time dimensions. 

Prints maximum of 10 columns per page. Number of 
rows printed is the same as the number of degrees 
of freedom. NRA is the dimensioned number of rows 
of the matrix and M and N are the number of rows 
and columns respectively of the matrix to be printed. 

MOUTC Used to print real and imaginary components of com- 
plex matrices. 

Uses object time dimensions. 

Prints maximum of 5 columns of real and imaginary 
data per page. NRA is the dimensioned number of 
rows of the matrix. The matrix to be printed is 
of order LxM. 

INVC Calculates the inverse of the complex impedance 
matrix to yield the complex mobility matrix. 

Uses object time dimensions. 

The original impedance matrix is destroyed. NRA 
is the dimensioned number of rows of the matrix 
and N is the order of the matrix. If the original 
matrix is singular IERR = -1 is returned to the 
main program, otherwise IERR = 0. 
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PROGRAM K2LRC 


INPUT 


Parameter 


Definition 


Card 


1 Columns 1 


2 


3 

4 


5 

6 


I Cl 


IC2 


IC3 

IC4 


IC5 

IC6 


Program Control 
= 1 Normal Program Execution 
= 9 Terminate Program 
1 Another Case Follows 

IC2 = 0 Applied Forces Only 
at Constrained Points 
IC2 > 0 Applied Forces at 
Non-Constraint Points in 
Addition to Applied Forces 
at Constraint Points 

Not Used 

>0 Print Simulated 
Measured Constrained Accel- 
eration and Simulated 
Measured Free Body Force 
Matrix 

= 0 Above Not Printed 
Not Used 


7-8 IC7 = 0 Use Mobility Data for 

Complete Frequency Spectrum 
from Tape IT1 Transmitted 
From Program KILRC 
>0 Use Mobility Data for 
Specific Frequencies as 
Indicated on Card, Maximum 
of 99 


9-80 HEADN Program Heading 
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Definition 


Card (s) 


Card (s ) 


Card 


Parameter 


FAMP Applied Force Vector 

One Force for Each Con- 
straint Format (8F10.0). 

Ten Columns Per Value, 

Eight Values Per Card, 
Maximum of 24 Applied Forces 

FAMP This card(s) Included Only 

IF IC2 > 0. Applied Forces 
at Non-Constraint Points. 
Format (8F10.0). Sum of 
the Number of Applied Forces 
at Constraint and Non- 
Constraint Points is Limited 
to 24 

Columns 1-5 PCT Uniformly Distributed (+) 

Random Percentage Error 
on Amplitude. Applied 
to Simulated Measured Con- 
strained Acceleration 


6-10 PCTB Constant Specified Per- 

centage Error. Applied to 
Simulated Measured Con- 
strained Acceleration 

11-15 PHE Uniformly Distributed (+) 

Phase Angle Error in 
Degrees. Applied to Simu- 
lated Measured Constrained 
Acceleration 

16-20 FPCT FPCT, FPCTB , FPHE are the 

same as PCT, PCTB, PHE 
Except Applied to Simulated 
Measured Forces 


21-25 FPCTB 


26-30 FPHE 

31-40 IZ Random Number Seed 
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Parameter 


Card 4 Columns 41-50 AMPL 


51-60 AMPR 


61-70 AMPF 


Card(s) 5 INDX 
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Definition 


Uniformly Distributed (+) 
Random Amplitude Error. 
(Simulating System Noise) 
Applied to Translational 
Acceleration 

Uniformly Distributed (+) 
Amplitude Error. (Simula- 
ting System Noise) Applied 
to Rotational Acceleration 

Uniformly Distributed (+) 
Random Amplitude Error. 
(Simulating System Noise 
Applied to Forces) 

This Card(s) Included Only 
IF IC7 > 0. Vector of 
Indices , in Ascending Order , 
Associated With Specific 
Frequencies , Indicating 
Which Corresponding Mobility 
Matrices Will be Used on 
Tape IT1 Transmitted From 
Program KILRC 

Format (8110) . Ten Columns 
Per Value, Eight Values 
Per Card 


defined. 



PROGRAM K2LRC - SUBROUTINES 


MOUTC Described in Program KILRC Subroutines . 

ERR2 Used to apply error to each element of a complex 

matrix. The errors include: a uniformly distributed 
random (+) percentage error on amplitude; a uniformly 
distributed (+) random phase angle; a constant 
specified percentage bias error; a uniformly distribu- 
ted (+) random amplitude error (simulating system 
» noiseT- Object time dimensions are used, with NRA 
the dimensioned number of rows. N1 and N2 are the 
number of rows and columns respectively of the 
matrix. 

GETRAN This subroutine is from the LRC Library of Subroutines. 
The subroutine computes uniformly distributed random 
real numbers between 0 and 1.0 (For Operation on the 
CDC 6400 and 6600 Computers) . 
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PROGRAM K3LRC 
INPUT, 


Parameter Definition 


Card 

1 Column 

1 . 

IC1 

Program Control 
= 1 Normal Program Execution 
= 9 Terminate Program 
^ 1 Another Case Follows 

* 

- ■ 

2 

IC2 

= 0 Standard Printed Output, 
Only Free Body Data 
= 1 Print Optional Output, 
Constrained Response Data 
= 2 Print Optional Output 
Force Matrix and Inverse 



3 / 

IC3 

Not Used 



4 

IC4 

IC4 ^ 0 Modify Response and 
Force Data on Tape IT1 
Transmitted From Program 
K2LRC 

IC4 - 0 No Modification 
of Above Data 



5 

IC5 1 




6 

IC6 

• Not Used 



7 

IC7 

4 




8 

IC8 

IC8 0 Print Complex Force 

and Response Data in Real 
and Imaginary Component 
Form 

IC8 = 0 Print Above Data 
in Amplitude and Phase 
Angle Form 



9-80 

HEADC 

Program Heading 
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Card 


Card 


Card 


Parameter 


Definition 


2 Column 1-80 ICHC Read only if IC4 ^ 0 Vector 

Indicating Which Response 
and Force Data to be 
Eliminated; One Element 
Per Channel 

ICHC = 0 No Modification 
ICHC = 1 Corresponding Data 
are Deleted 

Maximum of 80. Format 
(8011) . 


3 IW Read only if IC4 ^ 0 Vector 

Indicating Which Response 
Data to be Eliminated 
IW = 0 No Modification 
IW = 1 Corresponding Data 
are Deleted. One Element 
of IW for Each Force 
Eliminated. 

Maximum of 24 Format 
(2411) . 

4 Repeat Card 1 for options defined. 
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PROGRAM K3LRC 


(INPUT TAPE DESCRIPTION) 


This program converts constrained test data to free-body 
responses. The program accepts tape input in a form which 
would be generated from actual test data after processing. 
The input tape is unformatted and contains the following 
records of information. 


Record 

1 


2 


Parameters Description 

HEADT 80 characters (10 words) of heading 

NC Number of channels of data. Maximum 

of 80, 60 of which may be response 
data 

NF Number of channels containing force 

data (24 maximum) ; also the number 
of sets of data at each frequency. 
The remaining (NC - NF) channels 
contain response data (60 maximum) 


ICH Array (NC elements) , one element 

corresponding to each channel. 

0 indicates channel contains 
response data. 

2 indicates channel contains force 
data. 

1 added in the program indicates 
data to be ignored. 


FREQ Frequency at which NF records were 

taken. Identifies data which 
follows . 


3 W 

(Repeated 
NF times) 


Complex array (NC complex elements) , 
each element represents complex 
response or force data. There are 
as many W arrays as number of forces 
(NF) . Each W represents data taken 
under a single condition of forcing. 


Records 2 and 3 are repeated until FREQ = -1 to signal end of 
condition followed by next RECORD l. Final condition on 
tape is signalled by FREQ = -2. 
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PROGRAM K3LRC - SUBROUTINES 


MMPYC 


MOUTC 
INVC 
AMP HAS 


Performs multiplication of two complex matrices to 
yield a complex matrix as the result (C = A*B) . A 
is of order (LxM) . B is of order (MxN) and the 
resulting matrix C is of order (LxN) . Object time 
dimensions with NRA, NRB and NRC being the dimen- 
sioned rows in A, B and C, respectively. 

Defined in Program KlLRC Subroutines. 

Defined in Program KlLRC Subroutines . 

Converts the real and imaginary components of a 
complex matrix into amplitude and phase angle form. 
Object time dimensions with NRA being the dimensioned 
numbers of rows of the matrix. NROW and NFA are the 
actual number of rows and columns respectively of 
the matrix under consideration. 
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VARIED CONSTRAINTS PROGRAM 


2K1LRC , ZK2LRC , K3LRC 


These digital computer programs were designed to test and 
evaluate the concepts presented in the' varied constraints 
portion of the theoretical development. 5 The programs are 
similar to programs K1LRC and K2LRC. Both versions utilize 
program K3LRC to reduce the simulated constrained test 
data to the free body response of the structure. 

The same subroutines listed under program K1LRC and K2LRC 
apply. 
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PROGRAM ZK1LRC 


The input data for program ZK1LRC is the same as that for 
K1LRC with one exception. Card(s) 7 are repeated with the 
modified constraint characteristics. The number of sets 
of constraint data equals the number of constraints. 
Card(s) 7 are repeated one set for each constraint. 
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PROGRAM ZK2LRC 


The input data is the same as that for K2LRC with the 
exception of Card(s) 2. For program ZK2LRC , Card(s) 20 
are replaced with only one card having the following data: 

Parameter Definition 

Card 2 Column 1-10 F Applied Force 

11-20 NFS Coordinate Number of 

Constraint at Which 
Force is Applied 
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PROGRAM KlLRC 
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PROGRAM K2LRC 
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PROGRAM K2LRC 


PLOW CHART 
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PROGRAM K2LRC 


FLOW CHART 


SUBROUTINE 

MOUTC 

?RINT SIMULATED 
MEASURED 
CONSTRAINED 
ACCELERATION 


CALCULATE 
FORCES AT 
CONSTRAINTS 


SUBROUTINE 

ERR2 

APPLY ERROR 
TO FORCES 


/ HAS \ 
' NUMBER OF > 
FREQUENCIES 
. BEEN / 

\ reached/ 


/WRITE IT2 \ 
'acceleration' 

AND FORCE 
DATA FOR 
EACH 

^CONSTRAINT J 


SUBROUTINE 

MOUTC 

PRINT SIMULATE! 
MEASURED FREE 
BODY FORCE 
MATRIX 


f WRITE 
IT2 

FREQUENCY 
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PROGRAM K3LFC 
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PROGRAM K3LRC 
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PROGRAM K3LRC 
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PROGRAM 2K1LRC 


(VARIATION OF CONSTRAINTS) 
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PROGRAM ZKlLRC 


(VARIATION OP CONSTRAINTS) 
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PROGRAM ZK1LRC 


(VARIATION OF CONSTRAINTS) 
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PROGRAM ZK7LRC 


(VARIATION OP CONSTRAINTS) 
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PROGRAM ZK2LRC 
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o o o o o o o 


PROGRAM K1LRC < INPUT, OUTPUT, TAPE! ) 


NCON IS THE NUMBER OF CONSTRAINTS 
ONE INPUT CARO PER CONSTRAINT CONTAINING 

COORDINATE NUM6E R ,CONSTRA I NT ST I FFNESS , CONST R A I NT DAMPING 
NCOR IS THE INDICES OF THE STATIONS AT WHICH THE CONSTRAINTS ARE 
LOCATEO 

COMPLEX Z2, Y.Zi 

DIMENSION HE ADC I 9 ) , AM 1 60 . 61 I , G 1 60 I ,Z1 1 60.61 1 , GS I 60 ) , AM 60. 61) , 

1 NCOR (24) ,AKR(24) ,OAMP(24 ) ,GR (24 ) ,HZ 1 100 l . 

21 ROW (61) , I COL (61 ) , I RL AB ( 6 1 ) , I CL AB< 61 ) , Y 1 63, 6 1 ) , 1 2 t 60, 6 1 ) 

C NR A ON NEXT CARD MUST BE THE DIMENSIONED NUMBER OF ROWS OF MATRICES 

NR A *60 
I T 1 =1 

C READ FIRST CARD 

1 READ 1000,ICl,lC2,IC3,IC4,IC5,IC6,IC7,IC8,HEADC 

1000 FORMAT (811, 9A8) 

IF (IC1.EC.9) CALL EXIT 
PRINT 1005 

1005 FORMAT (54H1 K1LRC RAMAN AEROSPACE CORPORATION NOV 27, 1972 ) 

PRINT 1006, HE ADC 

1006 FORMAT (//25X.9A8//) 

1FIIC1 ,EU. 1) GO TO 100 

C READ FREE STRUCTURE DATA 

READ 1010, NO, GC 

1010 FORMAT (I10.7E10.0) 

C MASS MATRIX 

IFI.I.C 5 .£0. 1) GO TO 10 

1002 FORMAT (8E10.0) 

C DIAGONAL MASS 

DO 5 1*1, ND 
DO 5 J * 1 , NO 
5 AM ( I , J ) *0 • 

READ 1002, I A M ( 1 ,1 ) ,1*1, ND ) 

GO TO 20 

C FULL MASS MATRIX 

10 DO 15 1*1 ,ND 

15 READ 1002, (AMU, J),J*l,l) 

CALL SYM (AM, NO, NRA ) 

20 PRINT 1103 

1100 FORMAT (1HI/40X.11HMASS MATRIX//) 

CALL M0UT2 ( AM, ND,ND,NRA ) 

C R INPUT 

DO 25 1*1, ND 

25 READ 1C02, ( AR ( I , J ) , J* 1 , I ) 

CALL SYM ( AR , ND, NRA ) 

PRINT 120C 

1230 FORMAT ( IH1 /40X, 16HSTI FFNESS MATRIX//) 

CALL M0UT2 ( AR, ND , ND , NRA ) 

C VARIABLE DAMPING 

IF ( IC6 .CO. 0) GO TO 50 
READ 1002, I G < I ) ,I=1«NDI 
PRINT 1011, (G(I) ,1=1, NO) 
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-1011 FORMAT ( 1 OX . 1 3HDAMPI NG COEFS / ( 10X , 10F 10 . 3 I I 
00 30 I a 1|N0 
30 GSII) * SORT I C ( I-) ) 

DO 40 I > l.ND 

00 40 J > l.ND /- 

40 ailfj) * CMPLX (AMI , J) ,AK( I, J)*GS( I )*GS( JJ ) 

GO TO 100 

C CONSTANT STRUCTURAL DAMPING 

50 00 60 I « UNO 
DO 60 J - It NO 

60 Zl(ItJ) * CMPLX (AMI tJ) . GC*AM I , J ) ) i 

PRINT 10Z5t GC 

1025 FORMAT ( / /I OX f27H STRUCTURAL DAMPING COEFF. », F6.3I 
C CONSTRAINT DATA 

100 READ 1010. NCON 
00 180 I *1. NCON 

180 READ 1010. NCORI I ) . AKR 1 1 ) ,GRI 1 ) .DAMP 1 1 ) 

PRINT 1040, INCORII). AKR(I), GRID. 0AMPII),I-1, NCON) 

10 AO FORMAT I 1H1//10X.70HC00R0INATE SPRING RATE STR OAMPIN 

1G VI SCOUS DAMPING ///( 1 16, 1 P3E20.3 ) ) 

00 110 1 - l.NCON 
K > NCORI I) 

110 Z1IK.M * I1IM) ♦ CMPLX (AXRII), GRI 1 )*AKR( I ) ) 

FREQUENCY DATA 

READ 1010. NF 

READ 1002. (HZII I ,1-l.NF) 

WRITE tlTl) HEADC.NF.NO.IHZm.I-l.NF) 

WRITE (IT1I NCON, (NCOR(I) .AKRm.OAMPI I). GRID. I«l. NCON ) 

00 170 L»l.NF 
OMR*HZ I L 1*6.283185 
OMRS “OMR* OMR 
IFIICS.EQ.O) GO TO 140 
00 130 I *1 .NO 
00 130 J*i.ND 

130 Z2(I,J)*Z1(I » JI-OMRS*AMI I .0) 

GO TO 120 
140 00 135 I-l.NO 
00 135 J-l.NO 
135 Z2(I,J>-Z1(I,J> 

00 150 I-l.NO 

150 Z2U.n-Z2(I.I)-OMRS*AM(I,It 
120 00 160 1*1, NCON 

K*NCOR( I ) 

160 Z2(K.M*Z2(K,K> + CMPLX! 0.0 ,OMR*DAMPII) ) 

CALL INVC ( Z2.Y.N0 .NRA.NRA.IERR, I ROW, I COL ) 

IF i IERR.EQ.0 ) GO TO 151 
PRINT 1015 

1015 FORMAT (IH1/48H IMPEOANCE MATRIX (S SINGULAR JOB TERMINATED/ I 
CALL EXIT 

C PRINTED OLTPUT IS DISPLACEMENT OATA 
C 

151 PRINT 1017, HZIL) 
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1017 FORMAT! 1H1/A0X»A8HDI SPLACEMENT MOBILITY REAL , IMAG INARY 

1F10.2.6H HERT///1 

CALL MOUTC I Y,ND,ND,NRA, 0, 0 , I RL AB, ICL AO, 7 > 

190 WRITE II Til H 2 i L ) » ( I Y I I , J ) . 1 * 1 • NO ),J»1,N0> 

170 CONTINUE 
GO TO 1 
C 

END 
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SUBROUTINE SVM IA.N.NRA ) 

C FORMS SYMMETRIC MATRIX FROM LOWER TRIANGLE 

DIMENSION A I NR A« 1 I 
Nl-N-l 

00 13 l-I.NI 
1 1 -I ♦ X 

00 10 J«ll*N 
10 A(I,J)-A( J.l) 

RETURN 

ENO 
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SUBROUTINE MOUT2 (A ,M,N,NRA t 
DIMENSION 
1 A(NRA,1 t 

ID-MINOIN.IO) 

PRINT 1000, (1,1*1,10) 

1000 FORMAT l/SX, 10112) 

PRINT 100C 

DO 10 1*1 ,M 

10 PRINT 1001, I , I A ( I , J) , J* 1 ,10) 

1001 FORMAT (I5,5X,1P10E12.4) 

IF (IO-N) 20,50,50 

20 K-N/10-1 
00 40 1*1, K 
Nl*L* 10* 1 
N2*10* ( !♦ 1 ) 

ID-MI N0( N,N2 ) 

PRINT 1000, (1,1 *N1 , ID ) 

PRINT 1000 
DO 30 I *1 ,M 

30 PRINT 1001, I , ( A ( I , J ) , J* N1 , 1 D ) 
40 CONTINUE 

I F ( N2-N ) 60,50,50 
60 N2*N2*1 

PRINT 1000, 1 1 , I*N2 ,N) 

PRINT 1000 
DO 70 1*1, M 

70 PRINT 1001, I , (All , J) ,J-*2,N ) 
50 RETURN 
ENO 


93 



o o o o o o n 


subroutine moutc (a,l,m,nra,ir,ic, irlab, iclab, line » 

A IS COMPLEX ARRAY 1*16 ON IBM ) LXM 
IRLAB IS ARRAY OF INDICES FOR ROW (USED WHEN IR-I > 
ICLAB IS ARRAY OF INDICES FOR COL (USED WHEN IC«1> 
NRA IS DIMENSIONED NO OF ROWS IN A 
OUTPUT FORM RE.IM X.XXXE XX* X.XXXE XX 

LINE IS LINE NO ON PAGE OF FIRST OUTPUT LINE 

COMPLEX A (NRA, 1) 

DIMENSION IRLABI 1 ) * ICLABI II 
IC 1-1 

1C2-MIN0I5.M) 


10 

IL 1-1 




IL2-MIN0 (55 

-LINE «L ) 

15 

IF (IC 

.EQ.l) 

GO TO 20 


PRINT 

1000, 

(I,I-IC1,IC2 ) 


GO TO 

30 


20 

PRINT 

1000, 

(ICLAB (I) ,1— IC1.IC2 ) 

1000 

FORMAT 

(/ I23.AI2A/I 

30 

DO 50 

l-ILl, 

IL2 


IF IIR .EQ.l ) GO TO AO 

PRINT 1010* I . (A ( I ,J> .J-IC1.1C2) 

GO TO 50 

AO PRINT 1010 . IRLABII),tAII.JI,J-ICl,iC2 I 
1010 FORMAT I I6,AX,lP5(E12.3,lH,,Ell.3n 
50 CONTINUE 

IF (I.GE.L I GO TO 100 
IL1*IL2*1 

1L2<*MIN0( 1L2+5S.L ) 

PRINT 1015 
1015 FORMAT (1H1 I 
GO TO 15 

100 IF (IC2.GE.M ) GO TO 120 
IC1-IC2*1 

IC2-MIN0 ( IC2*5,M I 
PRINT 1015 
GO TO 10 
120 RETURN 
END 
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SUBROUTINE INVC ID ,A, N.NRA.NRB, IERR, IROW.ICOL ) 

C A - INVERSE Of 0 ORIGINAL MATRIX 0 IS DESTROYED 

COMPLEX A .0 « AMAX «E , PVT . AMULT 
OIMENSION AINRB.l) .DINRA.l) , I ROW III , I COL ( 1 ) 

I ERR *0 
DO 1 I -I, N 
DO 1 J-l.N 
1 AIl.JI-OM.J) 

H-N*l 

DO T I-l.N 
1R0WUI-I 

7 ICOL(I)-l 
DO 20 K-l.N 
AMAX- AIK.K) 

00 10 I-K.N 
00 10 J-K.N 

IEI CABS! AII.JII- CABS(AMAX) 110.9,9 
9 AMAX- All .Jl 
IC-1 
JC-J 

10 CONTINUE 
Kl-ICOLUI 

I COL I K) -I COL I IC) 

ICOLCIO-KI 

K I -I ROW I X ) f 

IRONIKI-IROHI JC» 

I ROW I JC ) -Kl 

IF( CABS! AMAX) ) 11,12*11 
12 1ERR — 1 

GO TO 100 

11 DO 14 J-l ,N 
E-AIK, J) 

A(K,JI«A(IC,J) 

14 AIIC.JI-E 
DO 15 I-l.N 
E-AI I,K) 

AII,K)*AII»JC) 

15 AII.JCI-E 
00 16 I-l.N 
IFI1-K) 18,17,18 

17 A(I,M)-(1.,0.) 

GO TO 16 

IB AII,M)«(0.,0.) 

16 CONTI NIC 
PVT -A (K,K ) 

DO 8 J-l.M 

8 A(K,JI-A(K,J) /PVT 
DO 19 I-l ,N 

IF I I-KI21 ,19,21 

21 AMULT-AII ,K) 

00 22 J-l.M 

22 A 1 1 , J ) -A I I »J) -AMULT* A ( K, J I 
19 CONTINUE 

00 20 I-l.N 
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20 A(l,Kt-A(t,*| 

00 25 1*1, N 
00 24 l-i,N 
IHIROWII >-024,23,24 

24 CONTINUE 

23 DO 25 J-l.N 

25 OIL, J)*A( I ,J> 

00 26 J*1,N 
00 28 L*1,N 

IFI ICOLf Jl-LI 26,29,28 

28 CONTINUE 

29 00 26 I *1 *N 

26 All, L) *01 I , J) 

100 RETURN 

END 
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PROGRAM K2LRC I INPUT .OUTPUT, TAPEI.TAPE2 1 
FREE BODY TEST METHOO SIMULATION 

USES MOBILITY OATA FROM ‘K1LRC* 


IC6 IS NUMBER OF DEGREES OF FREEDOM AT EACH COORDINATE 
USED IN SUBROUTINE ERR2 
COMPLEX Y , YA 

DIMENSION HEAONI 9), HEAD! 9 ) , I NOS C 2* I , INOX 1 100 I 

DIMENSION HZI100) .SPRI24) , FAMPI2A ) , DAMP! 24 I, ICHI 72), INDRI 20) 

DIMENSION Y(60,61),YA(60«61),GT(24),IRLAB(61),ICLAB(61) 

COMMON/SE T/LRCL,LRCN«LRCIX(2) 

TAPE UNIT ITI CONTAINS OATA FROM K1LRC PROGRAM 

TAPE UNIT IT2 USEO FOR OATA TRANSMITTAL TO K3LRC PROGRAM 


IT1-1 

IT2-2 

NRA-60 

FR«-1 

REMIND IT2 

C 

C READ FIRST CARO 

1 READ 1000,IC1,IC2,IC3,IC4,IC5,IC6,IC7, HE AON 
1000 FORMAT I6I1.I2.9A8 I 

IF I IC1.E0.9) CALL EXIT 
PRINT 1002 

1002 FORMAT (54H1 K2LRC KAMAN AEROSPACE CORPORATION NOV 27, 1972 I 
IF (IC1.EU.1) REMINO ITI 
15 READ (ITI) HEAO.NF.ND.IHZIIl.Ixl.NF) 

PRINT 1011, HE AON, HE AD, NO, (HZ(I),I B 1,NF) 

1011 FORMAT (1H1//25X.9A8// 

125X, 9AB /25X.I2.19H OEGREES OF FREEDOM /25X.23H FREOUENC 1ESI HZ I 
20N TAPE //(10X, 10F10, 2) ) 

C TAPE INPUT 

C CONSTRAINT SPRINGS, VI SCOUS DAMPERS, STRUCTURAL DAMPING 

REAO (ITI) NS • 1 1 NOS (I I ,SPR( 1 1 .DAMP! II.GTI I), 1*1, NS I 
C INPUT APPLIED FORCES AT CONSTRAINTS 

READ 130 , IFAMPII ),I«1, NS) 

C INPUT APPLIED FORCES AT NON-CONSTRAINT COROINATES 

IF I IC2.EQ.0) GO TO 7 
NS1-NS+1 
NS«IC2*NS 

READ 130, IFAMPII) ,I«NS1, NS I 
REAO 131, IINDSII ) ,1-NSl.NS I 
131 FORMAT (8110 ) 

130 FORMAT (8F10.0 I 
DO 11 l-NSl.NS 
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SPR ( 1 )=0. 

OAMP ( l I -0 . 

11 GTI I ) -O. 

7 PRINT 1003, UNOSI I) ,1-1, NS) 

PRINT 100A, I SPRII ) ,1 -1 ,NS> 

PRINT 1005, ( FAMPI I ) , I *1 ,NS ) 

PRINT 1001 , IDAMP(I) ,1-1 ,NS) 

PRINT 1006, (GT(I),I-1,NS ) 

1001 FORMAT I//10X,7HDAMPERS,10X,3I8F10.0)//) 

1003 FORMAT I / 10X.22HCONSTRA1NTS ANO FORCES // 10X , 9HSTAT ION ,5X, 

13(8110)//) 

100A .FORMAT I//10X, 6HSPRING , 1 1 X , 3 (6F10.0 ) // ) 

1005 format i // iox, unapplied force ,ax,3(8fio.oi// ) 

1008 FORMAT I / /10X , IBHSTRUCTUR AL DAMPING, 2X , 3 ( 8F10 .3 > // > 

C INPUT ERRORS, FREQUENCIES 

READ 1006, PCT ,PCTB , PHE iFPCT , FPCTB, FPHE, II , AMPL , AMPR , AMPF 

1006 FORMAT ( 6F5. 0 , I 1 0.3F10. 0) 

LRC I X ( 1 ) - I 

LRCN-1 

LRCL-1 

PRINT 1007, PCT, PCTB,PHE , 12»FPCT , FPCTB, FPHE, AMPL, AMPR, AMPF 

1007 FORMAT I/10X.16HMAX RAND ERROR «,F5.3,13H BIAS ERROR «,F5. 3,37m OF 

1 RESPONSE MAX RAND PHASE ERROR - ,F5.2,16H DEG. SEEO -» 110/ 

210X, 1 2HF0RCE ERRORS, AX .F5.3.13X ,F5.3, 37X.F5 .2// 10X, 3BHMAX RAND II 
3NEAR ACCELIREAL, I MAG) ERROR-E 12 .A/9X , 39HMAX RAND ANGULAR ACCELIREAL 
A, IMAGIERR0R-E12.A/17X.31HMAX RAND FORCEIREAL, IMAG ) ERROR-E 1 2. A// ) 

IF (IC7.GT.0) GO TO 20 
IC7-NF 

DO 10 1-1, NF 
10 !NDX(I)-I 

PRINT 1009. NF 

1009 FORMAT I/10X, AHALL ,13, IX, 16HFREQUENCI ES USED) 

GO TO 30 

20 READ 1OA0, 1 1 NO X ( I ) ,1-1, IC7) 

DO 12 1*1 ,IC7 
K* I NOX 1 1 ) 

12 Hill) -HI IK) 

30 PRINT 1010, (HZ I K) , K*1 , 1 C7 ) 

1010 FORMAT I /10X.16HFRE QUENC l ES USED// 1 10X , 10F12 .A ) I 
INFR-1 

DO 2 1*1, NO 

2 ICHII )-0 
NB-NOn 
NC«ND*NS 

00 3 I-NB.NC 

3 ICHI l ) *2 

WRITE IIT2) HEADN,NC«NS,( ICH(I ) , I * 1 , NC ) 

C START MAIN LOOP 

DO 500 L*1 »NF 

READ (ITl) FREQ, I I VII ,JI , I - 1 , NO ) , J- l , NO ) 

IFIL.NE.INDXIINFR)) GO TO 500 

C ELIMINATE COLUMNS AND CONVERT TO ACCL MOB 

OMR-FREQ *6.2832 
OMRS -OMR* OMR 
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DO 50 J*1,NS 
F=FAMP(J) 

K* I ND S ( J ) 

DO 50 I*1,ND 
Y(I,J)-V(I,K)*F 
50 YACl ,J)— V(| ,J)*0MRS 

C FORM CONSTRAINT OISP MATRIX ANO ADO ERROR 

CALL ERR2 ( YA , PCT ,PCTB , PHE . ND.NS , I X, NR A, AMPL , AMPR , 0 . » IC6 ) 

C 

IF I IC4.EQ.0t GO TO 60 
C 

PRINT 1020, H2IINFRI 

1020 FORMAT (1H1«20X,.71HSI MULATED MEASURED CONSTRAINED ACCELERATION 
l F - , F10 .2 t 3H HI/ t 

CALL MOUTC I YA ,ND,NS,NRA, 0, 0 , l RL AB, I CL AB, 5 ) 

ACCELERATIONS OUE TO FORCES ANO FORCES AT CONSTRAINTS 
WITH ERRORS ON TOTAL FORCES 
ACCELERATIONS OUE TO FORCES 

60 00 70 J-l.NS 
C FORCES AT CONSTRAINTS 

00 70 11*1, NS 
I-INDSlin 

70 Y(I1,J)*~ SPR (Il)*Y(I,J)-(GTIII)*SPR(II) ♦QNR*DAMP( II) l*Y( I.JI* 

1 10., 1.) 

DO 91 1*1 ,NS 
91 YII.M-YIl.lt+FAMPUl 

CALL ERR2 ( Y,FPCT ,F PCTB , FPHE , NS , NS . I X,NRA, 0. , 0 . , AMPF , 0 ) 

C 

IF ( IC4.EQ.0) GO TO 100 
C 

PRINT 1025, HZttNFR) 

1025 FORMAT ( 1H1 , 20X.60HSI MULATEO MEASURED FREE BODY FORCE MATRIX 
1 F - , F10.2.3H HI/ ) 

CALL MOUTC ( Y.NS.NS.NRA, 0, 0, l RLAB, ICLAB, 5 ) 

100 WRITE (IT2 I HZIINFR) 

C 

C ACCELERATION DATA ON TAPE 

9 DO 5 J-1,NS 

5 WRITE (IT2 ) (YA(I,J),!-1,N0),<YU,J>,!-1,NS ) 

1040 FORMAT (8110 ) 

502 INFR-INFR*1 

IF ( INFR.GT. IC7) GO TO 501 

500 CONTINUE 

501 WRITE ( I T2 I FR 
13 GO TO 1 

END 
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SUBROUTINE MOUTC ( A ,L, M.NRA , I R, I C, I RL AB, I CL AB, L INE > 

i 

A IS COMPLEX ARRAY 1*16 ON IBM ) LXM 
IRLAB IS ARRAY OF INDICES FOR ROM (USED MHEN |R*l t 
ICLAB IS ARRAY OF INOICES FOR COL (USED MHEN IC«1> 
NR A IS DIMENSIONED NO OF ROMS IN A 
OUTPUT FORM RE.IM X.XXXE XX. X.XXXE XX 

LINE IS LINE NO ON PAGE OF FIRST OUTPUT LINE 

COMPLEX A (NRA.li 
DIMENSION IRLAB(i) . ICLABI II 
ICI*1 

IC2«MIN0(5,M> 

10 I L 1*1 

IL2-MIN0 (55-LINE, L ) 

15 IF IIC.EQ.1) GO TO 20 

PRINT 1000, 1 1 ,I*(C1, IC2 I 
GO TO 30 

20 PRINT 1000, ( ICLAB (I) ,I*IC1,IC2 ) 

1000 FORMAT (/ 123,4124/1 
30 DO 50 I-IL1.IL2 

IF (IR .EQ.l ) GO TO 40 

PRINT 1010, I , (A II , J) ,J*1C1,1C2) 

GO TO 50 

40 PRINT 1010 , IRLAB(I),|A(I,JI,J*IC1,IC2 ) 

1010 FORMAT I I6,4X,1P5(E12.3,1H,,E11.3>) 

50 CONTINUE 

IF (I.GE.L I GO TO 100 
IL1-IL2»1 

IL2-MIN0I 1L2»55,L I 
PRINT 1015 
1015 FORMAT I 1H1 I 
GO TO 15 

100 IF (IC2.GE.M ) GO TO 120 
1C1*IC2+1 

IC2-MIN0 I IG2+5, M ) 

PRINT 1015 
GO TO 10 
120 RETURN 
END 
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SUBROUTINE ERR2 ( 24 ,PC T , PCTB , PHE , Ni , N2 . I X , NR, AMPL , AMPR . AMP , | C 6 ) 
OBJECT TIME DIMENSIONS 

EACH ELEMENT OF A COMPLEX MATRIX, A IS MODIFIED TO 
INCLUOE A SMALL PHASE ERROR, PNE ( DEG ) , A BIAS ERROR. 
PCTB (RATIO) ON AMPLITUDE, A UNIFORM RANDOM ERROR 
HAVING At/- MAXIMUM OF PCT (RATIO) ON AMPLITUDE, 

AND A 'UNIFORM RANDOM ERROR HAVING A ♦/- 
MAXIMUM ON AMPLITUDE 

THE PHASE ERROR IS ALSO RANDOMLY DISTRIBUTED 
NOTE NO SVMMETRI ZATION IS PERFORMED 

I ~ 6 IS THE NUMBER OF DEGREES OF FREEDOM OF EACH COORDINATE 

USES GETRANIUNIFORM DISTRIBUTION) 

LRC I X ( l ) * ARB I TRARI L Y SELECTED LOSITIVE INTEGER 
LRCN-1 FOR FIRST CALL TO GETRAN 

SET LRCN GREATER THAN 1 FOR SUBSEQUENT CALLS TO GETRAN 
LRCL»l FOR UNIFORM DISTRIBUTION 
LRC I X ( 2 ) » DUM1 , DUM2 , NOT USEO 
YFL-UNIFORMLY DISTRIBUTED RANDOM NUMBER 


COMPLEX 2A ( NR , i ) 

COMMON/SET/LRCL,LRCN,LRCI X(2) 

IF(PCT 1 120,100,120 
100 IF(PCTB) 120,110.120 
110 IF(PHE) 120,140,120 
140 IF ( AMPL ) 120,145,120 
145 l F ( AMPR ) 120,155,120 
155 I F ( AMP ) 120,135,123 
120 P=PHE /57. 296 
DO 130 J*1,N2 
K=(i:6tl)/2 
X1*IC6 

DO 130 1*1, Nl 

CALL GETRAN (LRC I X(1 ) , LRCN , LRCL , DUM 1 , YFL.DUM2) 

IFILRCN.EQ.l )LRCN*LRCNtl 
E«2.0*P*( YFL-0.5) 

AWAII.J) 

4I*(0.,-l.)*2A(I , J ) 

Ri* A3 S ( AR I 
R2* 4 B S ( A I ) 

R 3*R 1 *R2 

IF (R3.EQ.0.) GO TO 130 
Rl*Rl/R3 
R2*R2/R3 
Al*AR-E*Al 
A I *41 tSPAR 
AR.Al 

CALL 3ETRANILRCI X(U , LRCN ,LRCL , OUMi , YFL , DJM2 ) 

E*1.0t2.0*PCT*( YFL-0.5 I tPCTB 
A R * A R *c 
A I *4 I *fc 

CALL GETRAN(LKCl X ( l ) , LRCN , LRCL , DUMl , YFL , DJM2 ) 

C AMP ERROR DIST IN PROPORTION TO RATIOS OF REAL AND IMAG TO (REALt|MAG) 
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IF ( IC6.EQ. 0 I GO TO 170 
IF! l-K ) 176,175,185 
175 AMP-AMPL 
GO TO 170 
185 AMPaAMPft 

170 t«2.* CYFI-.5 I*AHP*R1 
AR«AR*E 

CAU GETRANILRCl XU ) , LRCN , LRCL , CXJM l , YFl , 0UM2 ) 

E«2.MVFL-.5 J*AMP*R2 

A1-A1*E 

160 ZAII.JI* CMPLX (AR.AI I 

IF I IC6.E0. 0 1 GO TO 130 
IF ( I .NE • K1 I GO TO 130 
K*IUIC6 
K1>KUIC6 
130 CONTINUE 
135 RETURN 
ENO 
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PROGRAM K3LRCI INPUT, OUTPUT, TAPE2) 

COMPLEX W,MF,WR,Y,HFINV,WF1 

01 MENS I ON H(80),HEA0C(9),HEA0T( 9) , ICHI80) , ICHCI 80), IWI24), 

1 WF(24,25),WR(60,24),Y(63,24) ,ICHNOI60),WFINV(24,25), ' 

2 IRON 125) ,!C0M25),HF1I24,25> 

C 

C TAPE UNIT I T1 CONTAINS DATA FROM F RE BOOT PROGRAM 

C 

ITl-2 

C READ FIRST CARO 

1 READ 1000, ICl,IC2,IC3,IC4,tC5,IC6,IC7, IC8, HEADC 
1000 FORMAT 1811, 9A8) 

IF I IC1.EQ.9) CALL EXIT 
PRINT 1005 

1005 FORMAT I54H1 K3LRC KAMAN AEROSPACE CORPORATION NOV 27 1972 ) 

IF I IC1.E0.1) REMIND IT1 

C READ FIRST TAPE RECORD 

READ IIT1) HEADT ,NC,NF,(ICH(I),I-1,NC) 

NIX*0 

NCA-NC 

NFA-NF 

C PRINT HEADINGS 

PRINT 101C, HEADC, HEAOT 

1010 FORMAT I//20X.9A8//20H TAPE HEADING , 9A8 ) 

PRINT 1015, NC,NF 

1015 FORMAT I /20X, 20HNUMBER OF CHANNELS » ,14, 10X , 18HNUMBER OF FORCES * 
1 ,14) 

IF IIC4.E0.0) GO TO 50 
C MODIFY TAPE OATA 

READ 1020 , ICHC ,1 M 
1020 FORMAT (8011) 

DO 10 I-l.NC 

IF IICHC(I).EO.O) GO TO 10 
NCA-NCA-1 

IF (ICHIl I.EQ.2) NFA-NFA-1 
ICHID-l 
10 CONTINUE 
NIX-NF-NFA 

IF INI X.EO.O) GO TO 30 
IX-O 

00 20 1-1,24 
IFIIMID.GT.O) IX-1XU 
20 CONTINUE 

IF ItX.EQ.NIX) GO TO 30 
PRINT 1025, NIX, IX 

1025 FORMA T I // 10X ,40H I NCORRECT NUMBER OF DATA SETS ELIMINATED , 15, 

1 9HREQUI REO, , I 5 , 10HDES I GNATED) 

CALL EXIT 

C PRINT MODIFICATIONS 

30 PRINT 103C.NCA ,NF A , (I, I M 1 1 I , I- 1 , NIX ) 

1030 FORMAT (20H MODIFIED TO ,20X, I4.28X, I4//20X, 34HDATA SETS EL 

1IMINATED - 1, KEPT - 0 // 7X, 12(15, 3H - , 1 2 ) /7X, 12 ( 1 5, 3H - ,12)) 

C CHANNEL INFORMATION * 

50 PRINT 1035, ( I , I CH ID , I -1 , NC ) 
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1035 FORMAT ( / /20X , 70HCHANNE L DEFINITIONS. - 0 RESP DATA, = 1 NOT U 

1SE0, 3 2 FORCE DATA // 7X.12U5.3H - , 12 1/ 7X, 121 15. 3H , ,121/ 

27X.12M5.3H = ,I2I/7X,12II5,3H = , 12 5/7X, 121 15, 3H 3 ,12)1 

C SET UP RESPONSE CHANNEL NOS 

NROW-O 
DO 60 I 3 1 • NC 

IF 1 1 CHI I l.NE.O) GO TO 60 
NROM 3 NR0W*1 
ICHNO(NROk) *1 
60 CONTINUE 

1FINR0W.GT.60) PRINT 1036 

1036 FORMATU0X.36H TOO MANY RESPONSE CHANNELS ) 

CALL EXIT 

310 READ IIT1) FREQ 

IF (FREO.GE.O) GO TO 90 
IF (FREQ. £0.-1) GO TO 1 
REWINO IT1 
GO TO 1 

C READ TAPE DATA 

90 NW-0 
ICOLO-O 
100 NW-NW*1 

IFINW.GT.NF) GO TO 200 
READ ( I T 1 I (Nil) .1-1, NC) 

IF INI X.EO.O) GO TO 120 
IFIIWINW) .NE.O) GO TO 100 
C DISTRIBUTE TAPE DATA 

120 1C0L0-IC0L0*1 
IROWF -0 

IROWR-O i 

DO 150 I-l.NC 
IF (ICHII)-l) 130,150,160 
130 IR0WR-IR0kR*l 

WRI IROWR. ICOLO) • W ( I ) 

GO TO 150 
160 IROWF -IROWF* l 

WF I IROWF .ICOLO)-W(I) 

150 CONTINUE 
GO TO 100 

C ALL INPUT MATRICES REAO-IN INVERT WF , 

200 IF I IC2.NE .2) GO TO 250 
C PRINT FORCES 

IF! ICB.NE. 0) GO TO 360 
PRINT 1051, FREQ 

1051 FORMAT I 1H1 , ///10X, 12H FREQUENCY - G16.6//10X, 3 1HF0RCE MATRIX 

1AMP, PHASE (DEG) //) 

DO 361 I-l.NFA 
DO 361 J-l.NFA 
361 WF1I1 , J ) > WF { 1 , J) 

CALL AMPHAS ( WF1 , NFA.NFA ,26 ) 

CALL MOUTC ( WF l , NF A ,NF A ,26 ,0 .0 , I CHNO, I CHNO, 7 , I C8 ) 

GO TO 250 

360 PRINT 1050, FREQ 

1050 FORMAT (1H1.///10X, 12H FREQUENCY 3 G16.6//10X, 31HF0RCE MATRIX 


104 



IRE AL • IMAGINARY //) 

320 CALL MOUTC C WF ,NF A • NF A ,24 ,0 ,0 , ICMNO, I CHNO, 7, ICB I 
250 CALL INVC I WF , WF I NV ,NF A, 24 ,2* , I ERR, I ROW , I COL ) 

IF ( IERR.EQ.OI GO TO 213 
PRINT 1040 

1040 FORMAT (///10X.45H FORCE MATRIX NOT INVERTIBLE ) 

GO TO 1 

210 CALL MMPYC I WR ,WF I NV, Y .NROW.NF A, NF A, 60. 24, 63 I 
C CALCULATIONS COMPLETE 

C PRINT OUTPUT 

IF (IC2.NE.2) GO TO 260 
IF ( IC8.NE.0 I GO TO 261 
CALL AMPHAS ( WF I NV ,NFA,NFA,24 I 
261 PRINT 1055, FREQ 

1055 FORMAT I 1H1///10X«12H FREQUENCY * G14 .4// 10X , 23HINVER SE OF FORCE 

1MATRI X //) 

CALL MOUTC (WFINV,NFA,NFA,24,0,0,ICHNO,ICHNO,7,IC8 > 

260 IF (IC2.E0.0) GO TO 330 

C PRINT CONSTRAINED RESPONSES 

I F I IC8.NE. 0 I GO TO 270 
PRINT 1061, FREQ 

1061 FORMAT (1H1.///12H FREQUENCY * G14.4// 10X , 4 5HC0NS TRAINED RESPONSE 

1 MATRIX AMP, PHASE (DEG) //) 

CALL AMPHAS ( WR ,NROM,NFA ,60 ) 

GO TO 280 

270 PRINT 1060, FREQ 

1060 FORMAT (1H1,///12H FREQUENCY * G14.4//10X.45HCONSTRAINED RESPONSE 
1 MATRIX REAL, IMAGINARY //) 

280 CALL MOUTC ( WR.NROW ,NFA,60, 1 ,0 , ICHNO, ICHNO, 7, 1C8 I 
330 IF I ICS .NE. 0 ) GO TO 300 
PRINT 1071 ,FREQ 

1071 FORMAT (1H1.///12H FREQUENCY = G14.4//10X.43HFREE BODY RESPONSE M 

1ATRIX AMP, PHASE (DEG) //) 

CALL AMPHAS I Y,NROW ,NFA ,60 ) 

GO TO 290 

300 PRINT 1070, FREQ 

1070 FORMAT (1H1.///12H FREQUENCY - G14.4// 10X , 43HFREE BODY RESPONSE M 

1ATRIX REAL, IMAGINARY //) 

290 CALL MOUTC I Y.NROW, NF A ,60 ,1 ,0 • ICHNO, ICHNO, 7 • ICB ) 

GO TO 310 
END 
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SUBROUTINE HMPYC I A * B , C , L , M, N, NRA, NRB, NRC I 

A t B f C COMPLEX MATRICES («16 ON IBM > 

C = A ♦ B A IS 11 H, H IS MXN, 

NRA.NRB.NRC DIMENSIONED ROWS IN A.B.C 

COMPLEX A INRAf 1) • B(NRU.t), CINRC.l) 

00 100 I-l.L 
00 100 J-I.N 
ClltJI -(0.0, 0.0) 

DO 100 K-l.M 

100 CII.J) - CII,J) ♦ All » K) * BIK.J) 

RETURN 

ENO 


C IS L X N 
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SUBROUTINE MOUTC ( A ,L , M*NRA, I R, I C , I RL AB, ICL AB, L INE > 

A IS COMPLEX ARRAY 1*16 ON IBM ) LXM 
IRLAB IS ARRAY OP INDICES FOR ROM (USED WHEN |R»l ) 
ICLAB IS ARRAY OF 1NOICES FOR COL (USED WHEN IC«1> 
NR A IS DIMENSIONED NO OF ROMS IN A 
OUTPUT FORM RE.IM X.XXXE XX* X.XXXE XX 

LINE IS LINE NO ON PACE OF FIRST OUTPUT LINE 

COMPLEX A (NRA, II 
DIMENSION IRLA81 1 1 * ICLABI 1) 

IC 1-1 

IC2-MIN0I5.M) 

10 ILI-1 

1L2-MIN0 (55-LINE, L I 
15 IF (IC.EQ.l) GO TO 20 

PRINT 1000, (I,I-IC1,IC2 I 
CO TO 30 

20 PRINT 1000, (ICLAB(I) .I-IC1.1C2 I 
1000 FORMAT (/ 123,41 24/ I 
30 DO 50 I-IL1.IL2 

IF HR .EQ.l I GO TO 40 

PRINT 1010, I • t A ( 1 ,J) ,J>1C1,IC2) 

CO TO 50 

40 PRINT 1010 , IRLAB(I) ,(A(I,JI,J«IC1, IC2 I 
1010 FORMAT ( I6,4X,1P5(E12.3,IH,,E11.3I) 

50 CONTINUE 

IF (I.GE.L I GO TO 100 
ILI-IL2M 

IL2-MIN0I IL2+55.L ) 

PRINT 1015 
1015 FORMAT (lHl I 
GO TO 15 

100 IF (IC2.GE.M I GO TO 120 
IC1-IC2+1 

IC2-MIN0 ( I C2+5, M ) 

PRINT 1015 
GO TO 10 
120 RETURN 
END 
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SUBROUTINE I NvC ID ,A , N, NRA , NRB, I ERR, I ROW , ICOL ) 

C A • INVERSE OF 0 ORIGINAL MATRIX 0 IS DESTROYED 

COMPLEX A ,D , AMAX ,E , PVT .AMULT 
DIMENSION A I NRB. l» .DINRA.l) , I ROW 1 1 ) , ICOL 1 1 ) 

IERR-0 
DO 1 I-l.N 
00 l J-l.N 
1 A( l,J)-DI I.J) 

M-N*l 

00 T I-l.N 
IR0WUI-! 

7 1C0LII)-I 
00 20 K-L.N 
AMAX- AIK.K) 

00 10 I-K.N 
DO 10 J-K.N 

IF ( CABSI AII.JII- CABS! AMAX) 110.9,9 
9 AMAX- All ,J) 

IC-1 

JC-J 

10 CONTINUE 
Kl -ICOL I K ) 

ICOL IK) -I COL I IC) 

ICOL I IC I -Kl 
Kl -IROW(K) 

1R0MIK) -1 ROW! JC) 

IRONI JCl-KI 

IF I CABSI AMAX) I 11,12,11 
12 IERR — 1 

GO TO 100 

11 DO 14 J-l.N 
E-AIK, J) 

AIK.JI-AI IC.J) 

14 AIIC.JI-E 
00 15 I-l.N 
E-AII.K) 

A ( I , K ) -A 1 1 , JC ) 

15 AII.JCl-E 
00 16 1-1, N 
IFII-K) 18,17,18 

17 A(I,M)-(1.,0.) 

GO TO 16 

18 AII,M)-(0.,0.) 

16 CONTINUE 
PVT-AIK.K) 

DO 8 J-l • M 

8 AIK, J) -A IK, J) /PVT 
00 19 1-1 ,N 

IF 1 1 — K ) 2 1 ,19,21 

21 AMULT-AII ,K) 

DO 22 J-l.M 

22 All, J) -A (I , J) -AMULT * AIK, J ) 

19 CONTINUE 

00 20 l-l.N 
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20 Att,K)*AII,M) 

00 25 1*1, N 
00 25 L«1,N 
IFUROWI l J-U25,23,25 
25 CONTINUE 
23 DO 25 J*1 , N 

25 OIL, JI»Al I ,J) 

00 26 J*1 ,N 
00 28 L*1,N 
IFIICOLIJI-U 28,29,28 

28 CONTINUE 

29 00 26 1*1 ,N 

26 A(I,LI*0(I«Jt 
100 RETURN 

END 
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SUBROUTINE AMPHAS ( V > NROW * NF A * NR A > 

C CONVERTS REAL ANO I NAG I NARV COMPONENTS TO AMPlITUOE AND PHASE ANGLE 

COMPLEX V 

DIMENSION VtNRA.l ) 

DO 1 I a l » NROW 
00 1 J = 1 * NF A 
Al-VI I .-1. ) 

A R * Y ( I i J ) 

I F I AI.EQ.O.O.AND.AR.EQ.O.O I GO TO 2 
AM- SORT! AR*AR+AI*A1 I 
Al- ATAN2 (Al t AR 1*57.2956 
GO TO 1 
2 AM-O. 

1 Y I I , J I - CMPLXIAM.Al ) 

RETURN 

ENO 
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PROGRAM ZKlLRC (INPUT. OUTPUT, TAPE 1.TAPF2) 

VARIATION OF CONSTRAINT CHARACTERISTICS 

NCON IS THE NUMBER OF CONSTRAINTS 

CNE INPUT CARO PER CONSTRAINT CONTAINING 

COORDINATE NUMBER, CONSTRAINT ST IFFNESS.CONSTRA I NT DAMPlNo 

NCOR IS THE INUiCtS OF THE STATIONS AT WHICH THE CONSTRAINT* AKt 

LOCATED 


COMPLEX 22,Y,Z 4 ,Z2A,Z3 

DIMENSION HEAOCI 9) , AMI 60.61I.GI 60) ,21(60,61) ,GSI60) , AM 60 ,6 1 ) , 

1 NCOR (24, 24) «AKR(24, 241 .DAMP! 24, 2 A) ,GR(24,24) ,H2U00) , 

2 I ROM! 6 1 ) , I CUL (61), IRL ABI 6 1 ) . ICL AB (61 ) , YI60 ,6 1 ) , 2 2 1 60 , 61 1 ,2i ( 2 4) 
DIMENSION 22A( 60,61 ) 

C NR A ON NEXT CARC MUST BE THE DIMENSIONED NUMBER OF ROWS OF MATKlv.ES 

NR A *60 
IT 1*1 

C READ FIRST CARO 

1 READ 1000, 1C1, 1C2, 1C3.1C4, ICS, IC6,ICT,ICB,HEADC 
1000 FORMAT (811,9Ae) 

IF (ICl.Ew.9l CALL EXIT 
PRINT 1005 

1005 FORMAT (54H1 K1LRC KAMAN AEROSPACE CORPORATION DEC 14, 1972 ) 

PRINT 1006, HEACC 

1006 FORMAT I//25X.9A8// ) 

IF( IC1 .EU. 1) GO TO 100 

C READ FREE STRUCTURE DATA 

REAC 1010, NO, GC 
1010 FORMAT ( 1 10 , 7E 10*0 ) 

C MASS MATRIX 

IF ( ICS * EQ, 1) GO TO 10 
1002 FORMAT (8E10.0) 

C 01 AGONAL MASS 

DO 5 1*1, NU 
CO 5 J* 1 , NO 
5 AMI I,JI*0. 

REAC 1002, (AMI I, I 1,1*1, NO) 

GO TO 20 

C FULL MASS MATRIX 

10 DO 15 1*1, NO 

15 REAC 1U02,(AM(1,J) , J* 1 , I ) 

CALL STM ( AM, NO ,NRA ) 

20 PRINT 1100 

1100 FORMAT ( 1H1/40X, 11HMASS MATRIX//) 

CALL MUUT2 ( AM ,NC, NO , NR A ) 

C K INPUT 

DO 25 1*1, NO 

25 REAC 1002, ( AM I , J ) , J = 1 , 1 > 

CALL SYM ( AK.NU.NKA I 
PRINT 1200 

1200 FORMAT 4 1H1/4UX , 16HST IFFNE SS MATRIX//) 

CALL M0LT2 ( AK , ND, ND.NRA ) 
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VARIABLE DAMPING 
IF IIC6 .EO. 0) GO TO 50 
READ 1002, (GIH, 1*1, NO) 

PRINT 1011, IGU) ,1-1, ND) 

1011 FORMAT ( 1 CX , 1 3HDAMP I NG COEFS / 1 10X, 10F10. 3 I ) 

DO 30 I - 1 , NO 
30 GSm - SORT (Gill) 

00 AO I - 1 , NO 
DO 40 J » 1 ,ND 

40 21(I,J) » CMPLX I AK( I , J) ,AK( I , J) *GSI I )*GS( J I J 
GO TO 100 

C CONSTANT STRUCTURAL DAMPING 

50 DO 60 I - 1,N0 
DO 60 J - 1 ,ND 

60 21(1, J» - CMPLX ( AX ( I , J ) , GC+AK ( I , J ) ) 

PRINT 1025, GC 

1025 FORMAT I //10X.27HSTRUCTURAL DAMPING COEFF. », F6.3) 

C CONSTRAINT DATA 

100 PRINT 1050 

1050 FORMAT I 1H1 ///10X,21HCONSTRA! NT VARIATIONS///) 

READ 1010, NCON 
DO 125 J»1,NC0N 
DO 125 I -1, NCON 

125 READ 1010, INCOR 1 1 , J) , AKRI I , J) ,GR( I , J ) , DAMP! I, J ) ) 

00 110 J-l.NCON 

PRINT 1040, INCOR 1 1 • J) ,AKR 1 1 , J) ,GRI I , J) , DAMP I 1 , J ) , I* 1, NCON ) 

110 PRINT 1060 
1060 FORMAT (////) 

C FREQUENCY DATA 

READ 1010, NF 

READ 1002, <H2(l),l»l,NFI 

WRITE IIT1) HE ADC , NF ,N0, ( HZ ( I ) , I"1 , NF ) , NCON 

DO 170 L-l.NF 

WRtTEIITl) Hllll 

0MR-H2IL >*6,2831 85 

OMRS-OMRPOMR 

IFIIC5.E0.0) GO TO 140 

DO 130 I-l.NO 

DO 130 J-l.NO 

130 2211 .JI-Z1II ,J)-OMRS*AM(I ,J| 

GO TO 120 
140 DO 135 1-1, ND 
00 135 J«1,N0 
135 22(1. J)*Z1(I ,J) 

DO 150 1-1, NO 

150 22( 1 , 1 1-22(1 , l)-OMRS*AM(I , I) 

120 UO 170 1 1 «1, NCON 

WRITE 1 1 T1 1 (NCORII ,11), A KR (I, II) , DAMP I I, II), GR (1,11), 1*1, NCON ) 

DO 160 I-l.NCON 
K « NCOR I 1 , 1 1 ) 

23(1 >-Z2IK,K> 

160 22(K,KI - 22 ( K ,K ) ♦ CMPLX ( AKR( 1 , 1 1 ) • OMR*DAMP( I , II ) ♦ GR 11,(11* 
1AKRI I , 1 1 ) ) 

1040 FORMAT ( //10X ,70HCOORO I NATE SPRING RATE STR DAMPIN 


112 



1G VISCOUS DAMPING /// 1 1 16 , 1P3E20. 3 > /// » 

C 

C 

DO 152 I-l.ND 
DO 152 J*1 ,ND 
152 Z2A ( 1 • J) >£2( I ( J) 

CALL IN VC I Z2A,Y,ND,NRA,NRA,IERR, IROW, ICOL ) 

IF I IERR.EQ.0 ) GO TO 151 
PRINT 1015 

1015 FORMAT I1H1/48H IMPEDANCE MATRIX IS SINGULAR 
CALL EXIT 

C PRINTED OUTPUT IS DISPLACEMENT DATA 
C 

151 PRINT 101T, HZIL) 

1017 FORMA T 1 1H 1/A0X,68H0 1 SPLACEMENT MOBILITY REAL 

1F10.2.6H HERTZ// I 

CALL MOUTC I Y.NOtND.NRA* 0, 0, I Rl AB, ICL AB, 7 > 

00 162 I-l.NCON 
K ■ NCORII.II) 

162 Z2(K,K)-Z3III 

190 WRITE C I T 1 I I ( Yt I ,J) , I-l.ND >,J«1,ND) 

170 CONTINUE 
GO TO 1 
C 

END 


JOB TERMINATED/ I 


IMAGINARY FREQ- 


113 



SUBROUTINE SYM IA.N.NRA ) 

C FORMS SYMMETRIC MATRIX FROM LOWER TRIANGLE 

DIMENSION AINRA.l I 
Nl-N-1 

DO 13 1*1, Nl 
1 1 -I ♦ l 

00 10 J*il,N 
10 A(I«JI*A(JtIi 
RETURN 
END 
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SUBROUTINE MOUT2 (A.M.N.NRA I 
DIMENSION 
1 A (NR A , 1 ) 

IO«NINO(N,10> 

PRINT 1000. (1,1*1.10) 

1000 FORMAT (/SX. 10112) 

PRINT 100C 

DO 10 1*1. M 

10 PRINT 1001* I , (A ( I , J) , J=1 , ID) 

1001 FORMAT I I 5.5X.1P10E12.0) 

IF (ID-N) 20.50.90 

20 K-N/10-1 
00 00 L*1»K 
Nl*l*10^l 
N2*10* I L* 1 ) 

IO-MINOI N.N2 > 

PRINT 1000. tl.l-Nl.ID) 

PRINT 1000 
DO 30 1*1 .M 

30 PRINT 1001. I ,(A(I,J). J*N1,ID) 
00 CONTINUE 

I F ( N2-N I 60,50.50 
60 N2-N2«l 

PRINT 1000, I l , I*N2 ,N) 

PRINT 1000 
DO 70 1*1, M 

70 PRINT 1001. I ,(A(I ,J) .J-N2.N I 

50 RETURN 
ENO 
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SUBROUTINE MOUTC 1 4 ,L .M.NRA, I R, I C, I RL AB, ICL AB, L INE ) 

A IS COMPLEX ARRAY (*I6 ON IBM ) LXM 
IRLAB IS ARRAY OF INDICES FOR ROW (USED WHEN IA«1 ) 
ICLAB IS ARRAY OF INDICES FOR COL (USED WHEN IC*1> 
NR A IS DIMENSIONED NO OF ROWS IN A 
OUTPUT FORM RE • IM X.XXXE XX. X.XXXE XX 

LINE IS LINE NO ON PAGE OF FIRST OUTPUT LINE 

COMPLEX A (NR A .11 
DIMENSION IRLABI 1) . ICLABI 1) 

ICl-i 

IC2-MINOI 5,M) 

10 I L 1-1 

IL2-MIN0 (55-LINE *L I 
15 IF ( IC.EQ.l ) GO TO 20 

PRINT 1000* (I.I-IC1.IC2 I 
GO TO 30 i 

20 PRINT 1000* ( ICLAB ( I ) .I-IC1.IC2 I 
1000 FORMAT (/ 123.4124/) 

30 DO 50 I-IL1.IL2 

IF HR .EQ.l ) GO TO 40 

PRINT 1010* I .(All . J) *J«IC1*IC2) 

GO TO SO 

40 PRINT 1010 . IRLAB (I).(A(I.J).J*IC1.IC2 ) 

1010 FORMAT I I6.4X, IPS (£12. 3.1H, *£11.31) 

SO CONTINUE 

IF (I.GE.L I GO TO 100 
IL1*IL2*1 

IL2-MIN0I IL2+SS.L I 
PRINT 1015 
1015 FORMAT (1H1 I 
GO TO 15 

100 IF ( IC2.GE.M ) GO TO 120 
IC1»IC2*1 

IC2-MIN0 (IC2+5.M ) 

PRINT 1015 
GO TO 10 
120 RETURN 
END 
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SUBROUTINE I NVC ID ,A , N, NR A , NRB, I ERR, I ROW, ICOL ) 

C A * INVERSE OF 0 ORIGINAL MATRIX 0 IS OESTROVEO 

COMPLEX A <0 • AMAX ,E , PVT , AMULT 
DIMENSION AINR8.ll , 0 1 NR A , 1 1 ,1 ROW ( 1 9 , l C QL ( 1 1 
I ERR *0 
DO 1 I-l.N 

00 1 J-l.N 

1 All,J)«Dtl,J) 

M»N*1 

DO 7 I-l.N 

1 row ( i i«i 

7 1C0L I I )*l 
DO 20 K«1 ,N 
AMAX- A I K «K) 

DO 10 I-K.N 
00 10 J-K.N 

IF! CABS! A 1 1 , J) ) - CABS I AMAX) 110,9,9 
9 AMAX- All ,J) 

IC-l 

JC-J 

10 CONTINUE 
KI-ICOLIK) 

I COL ( K) * I COL I IC) 

ICOLt ICI-Kl 
KI-IROWIK) 

IRQW(K) -I ROW! JC) 

1R0WI JC)-K1 

IF I CABSI AMAX) ) 11,12,11 
12 IERR— 1 

GO TO 100 

11 DO 14 J-l.N 
E-AIK, J) 

AIK.JI-AI IC.J) 

14 AIIC.JI-E 
00 IS I-l.N 
E-AII.K) 

A(I,K)-AII,JC) 

15 AII,JC)-E 
DO 16 I-l.N 
IFII-K) 18, 17, IB 

17 A(I,M)-{1.,0.) 

GO TO 16 

18 AII,M)»|0.,0.) 

16 CONTINUE 
PVT-A (K,K ) 

DO 8 J-l • M 

8 AIK, J) -A IK, J) /PVT 
DO 19 1-1 »N 

IF ( I-K) 21 ,19,21 

21 AMULT-AII ,K) 

DO 22 J-l.M 

22 All, J) -A II , J) -AMULT*A IK , J ) 

19 CONTINUE 

00 20 I-l.N 
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20 Ail , K ) >A ( I , H ) 

00 25 I«1,N 
00 2A L B 1 »N 
IFIIROWII >-D2*,23,2* 

2 A CONTINUE 
23 00 25 J-l.N 

25 OIL, J»«AI I , J) 

00 26 J-t.N 
00 28 1*1, N 
IFIICOUJI-L) 28,29.28 

28 CONTINUE 

29 00 26 (-1.N 

26 A(I,U»0(t,JI 
100 RETURN 

ENO 



no oooooono nnonoooooo 


PROGRAM 2 K2LRC ( INPUT, OUTPUT , T AP E 1 .TAP E2 ) 

FREE dUUY TEST METHOD SIMULATION 

VARIATION CJf CONSTRAINT CHARACTERISTICS 

USES MOBILITY DATA FROM ‘ZKILRC* 


IC6 IS NUMBER OF DEGREES OF FREEDOM AT 
US EO IN SUBROUTINE ERR 2 
COMPLEX Y.YA.Yl 

01 MENS ION HEAONI 91, HEAD! 9),INDS(24I, 
01 MENS ION HZI100I, I CHI 841, INDR I 20) 
DIMENSION Y(6U,61),YA(60,61),IRLAB(6I) 
DIMENSION SPR(24,24I,DAMP<24,24),GT(24 
COMMON /SET /LRCL , LRCN, LRC IXI2) 


EACH COORDINATE 


INDXI 1001 

, ICLABI61) 

,24) , Yl (60,61) 


TAPE 

UN IT 

IT 1 

CONTAINS 

OATA 

TAPE 

UNIT 

IT2 

USED FOR 

OATA 


FROM K1LRC PROGRAM 
TRANSMITTAL TO K3LRC PROGRAM 


IT 1 *1 
IT2-2 
NRA-60 
FR = -1 

REMIND 1T2 


RE AC FIRST CARC 

1 READ 100U, ICl, IC2, IC3, IC4, ICf,IC6,lC7, HEADN 
1000 FORMAT (611,12, SAB I 

IF (ICi.Ew.9l CALL EXIT 
PRINT 1002 

1002 FORMAT (54H1 K2LRC KAMAN AEROSPACE CORPORATION OEC 20, 1972 ) 

IF ( ICl.Ew.ll REWIND IT1 
15 RE AC ( I T 1 ) HEAU.NF.ND, (HZ( II , I = 1 , NF ) ,NS 

PRINT lull, FEAON, HEAC, NO, (HZ( II ,I>1,NF) 

1011 FORMAT ( 1H1//25X,9A8// 

125*, 9AB /25a ,12, 19H DEGREES OF FREEDOM /25X.23H FREUUENCI E S( H2 I 
2 ON TAPE //(10X, 10F10.2II 
C INPUT APPLIEC FORCE AT CONSTRAINT 

REAC 140 , F.NFS 
140 FORMAT (FlU.O, l 10 I 
C INPUT ERRORS, FREQUENCIES 

REAC 1 0u6 » PCT,PCTB,PHE,FPCT,FPCTB,FPHE,I7,AMPL,AMPK,AMPF 

1006 FORMAT ( 6f 5. 0 , 1 1 0, 3F 10. 0 I 
LRC IX ( 11 = 12*2*1 

LRCN=1 
LR C L= 1 

PRINT 10w7, FLT, PCTB ,PHE , IZ, FPCT , F PC T B , FPHE ,AMPL ,AMPK,AMPF 

1007 FORMAT </10A,16FMAX RAND ERROR =»F5.3,13H BIAS ERROR =,F5.j,j7H at 

l RESPONSE MAX RAND PHASE ERROR = ,F5.2,16H DEG. SEEu = , 110/ 
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210X, 12HF0RCE ERRORS , AX ,F 5 . 3 , 1 3X ,F5. 3 . 37X , F5 . 2// 10X , 38HM Ax RAND U 
3NEAR ACCELIREAL, l HAG) E RROR-E 12 . A/9X, 39HMAX RAND ANGULAR ACCEL I REAL 
A,IMAG)ERRGR*E12.A/17X.31HMAX RAND FORCE (REAL* IMAG > ERROR-E 12. A/ / ) 

IF ( IC7.GT.OJ GO TO 20 
IC7-NF 

00 10 I-l.NF 
10 1N0K(IJ-I 

PRINT 1009, NF 

1009 FORMAT I/10X* AHALL ,I3,lX, 16HFRE0UENC I ES USED) 

GO TO 30 

20 READ 10AJ, I I NDX (I) , I - 1 , IC7 I 
DO 12 I-1.IC7 

K-INOXU i 
12 H2II)«HZ(KI 

30 PRINT 101C, (HZ(K) ,K>1 ,IC7 ) 

1010 FORMAT (/10X.16HFRE0UENCIES USED// ( 10X * 10F12.AJ) 

INFR-l 

DO 2 I -1, NO 

2 ICHUt-0 
NB»ND*1 
NC-ND*NS 

DO 3 I-NB.NC 

3 ICHli i-2 

WRITE (IT2) HEADN.NC ,NS , I 1 CHI l),l-l*NC ) 

C START MAIN LOOP 

00 500 L-l.NF 

IFIL.NE.INOX(INFR) ) GO TO 300 
C TAPE INPUT 

C CONSTRAINT SPRI NGS , VI SCOUS OANPERS. STRUCTURAL DAMPING 

READ IIT1 ) FREQ 
OMR-FREQ *6.2032 
OMRS-ONR*OMR 
DO 21 LL-l.NS 
PRINT 1016. FREQ 

1016 FORMAT (///39H CONSTRAINT CHARACTERISTICS FOR FREQ-. F7. 2. 3H HZ/) 
READ (IT1 I 1 1 NO S ( I I.SPRII.LL) .DAMP! I.LLI.GTI I.LL), I-l.NS ) 
IF(L.NE.l) GO TO 7 
PRINT 1003, (INDSII ), I-l.NS) 

PRINT 100A « I SPRI I.LL). I-l.NS I 
PRINT 1001 , I DAMP (I.LL). I-l.NS) 

7 READ 1171) I IY1 II ,J) .1-1 .NO) .J-l.ND ) 

C ELIMINATE COLUMNS ANO CONVERT TO ACCL MOB 

DO 21 I-l.ND 
YU,LL)»Y11I,NFS)*F 

21 YAd.LLI— Y(I ,LL)*OMRS 

1001 FORMAT I //10X.7HDAMPERS. 10X ,3 (8F10.0)// ) 

1003 FORMAT (//10X.9H STATION ,5X ,3 1 8 1 10 ) / / ) 

100A FORMAT (//10X, 6HSPRING , 1 IX , 3 C 8F10.0 ) // I 
PRINT 1005, F, NFS 

1005 FORMAT ( / /10X .13HAPPLI EO FORCE , AX ,F 10 . 1, 3X, lOHCOORO INATE, 1 3/ / ) 

C FORM CONSTRAINT DISP MATRIX AND ADD ERROR 

CALL ERR2 ( YA , PCT .PCTB.PHE , NO. NS. IX, NRA, AMPL • AMPR.O. . IC6 I 

IF UCA.EQ.O) GO TO 60 
PRINT 1020, HZIINFR) 
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non 


1020 FORMAT I IH1 ,20X, 71HSI MULATED MEASURED CONSTRAINED ACCELERATION 
1 F «■ , F 10 • 2 • 3H HZ/I 

CALL MOUTC ( YA , NO ,NS , NR A, 0, 0 , I RL AB, ICL AB, 5 > 

ACCELERATIONS OUE TO FORCES AND FORCES AT CONSTRAINTS 
Ml TH ERRORS ON TOTAL FORCES 
ACCELERATIONS OUE TO FORCES 

60 00 TO J -1 ,NS 
C FORCES AT CONSTRAINTS 

00 70 11*1, NS 
i I ” I NO SI 1 I ) 

‘ 70 Y( II, J| — SPRIIItJI *Y ( I , J ) - (GT ( I I , J I *SPR I 1 1, J ) ♦ OMR* DAMP I I I , J I )* 
lYtI,J}*<0.,l.l 
00 91 J-l.NS 
DO 91 I-l.NS 

IFI NFS.NE.INDSIII ) GO TO 91 
Yl l » J I *YI I , J I *F 
91 CONTINUE 

CALL ERR2 ( Y, FPC 7 ,FPC 78, FPHE, NS , NS , IX, NRA, 0., 0., AMPF, 0 I 

IF (IC6.EQ.0I GO TO 100 
PRINT 1025, HZUNFRI 

1025 FORMAT I lHl , 20X, 60HSI MULATED MEASURED FREE BODY FORCE MATRIX 
1 F * , F10.2.3H HZ/ I 

CALL MOUTC I Y.NS.NS.NRA, 0, 0 , I RL AB, I CL AB, 5 > 

100 WRITE (IT2 ) HZUNFRI 

C 

C ACCELERATION DATA ON TAPE 

9 DO 5 J«1,NS 

5 WRITE (1T2 I (YAII, Jl, 1*1, N0>, (Y(I, J), I-l.NS I 
1060 FORMAT 18110 I 
502 !NFft-!NFR«l 

IF I INFR.GT.IC71 GO TO 501 

500 CONTINUE 

501 WRITE IIT2 ) FR 
13 GO TO 1 

END . 
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o o o o o o o 


SUBROUTINE MOUTC I A ,L,M,NRA, IR, IC, IRLAB, ICLAB, L INE I 

A IS COMPLEX ARRAY 1*16 ON IBM ) LXM 
IRLAB IS ARRAY OP INDICES POR ROW (USED WHEN IR-1 I 
ICLAB IS ARRAY OF INDICES FOR COL (USED WHEN IC»1> 
NR A IS DIMENSIONED NO OP ROWS IN A 
OUTPUT FORf< RE.IM X.XXXE XX, X.XXXE XX 

LINE IS LINE NO ON PACE OF FIRST OUTPUT LINE 

COMPLEX A (NRA.l) 

DIMENSION I RLAB (II, ICLAB ( 1 I . 

ici-i , . , , 

IC2-MIN0C 5.M) 

10 ILl-l 

IL2*MIN0 (55-LINE, L ,1 
15 IP (IC.EO.ll GO TO 20 

PRINT 1000, (I , I • IC 1 • I C2 I 
CO TO 30 

20 PRINT 1000, (ICLAB(I) , I > I Cl , I C2 I , 

1000 FORMAT (/ 123,4124/1 
30 00 50 I-IL1.IL2 

IP UR .EQ.l I GO TO 40 

PRINT 1010, I , (A ( I , J) , J«IC1,IC2) 

GO TO 50 

40 PRINT 1010 , IRLAB(I) , ( A ( I , J ) , J* I C 1 , IC2 ) 

1010 FORMAT ( 16,4X,1P5(E12,3,IH,,EU.3I) 

50 CONTINUE 

IP (l.GE.L I GO TO 100 
IL1«IL2*1 

IL2“MIN0( IL2*55,L I 
PRINT 1015 
1015 FORMAT (1H1 ) 

GO TO 15 

100 IF ( IC2.GE.M ) GO TO 120 
IC1-IC2*1 

IC2-MIN0 IIC2+5.N I 
PRINT 1015 
GO TO 10 
120 RETURN 
ENO 
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o n n n n n n n n n ri o n n n n r> n n n 


SUBROUTINE' ERR2 I l A ) PC T , PCTB , PHE , N1 , N2 . I X , NR , AMPL , AMPR, AMP, [ C 6 ) 
OHJECT TIME DIMENSIONS 

EACH ELEMENT OF A COMPLEX MATRIX, A IS MODIFIED TO 
INCLUDE A SMALL PHASE ERROR, PNE IDEG), A BIAS ERROR, 
PCTB (RATIO) ON AMPLITUDE, A UNIFORM RANDOM ERROR 
HAVING A ♦/- MAXIMUM OF PCT (RATIO). ON AMPLITUDE, 

AND A UNIFORM RANDOM ERROR HAVING A ♦/- 
MAXIMUM ON AMPLITUDE 

THE PHASE ERROR IS ALSO RANDOMLY DISTRIBUTED 
NOTE NO SYMMETRl ZATION IS PERFORMED 

I C6 IS THE number of DEGREES of FREEDOM OF EACH COORDINATE 

USES GETRANI UNIFORM DISTRIBUTION) 

LRC I X 1 1 1 » ARB I TRARI L Y SELECTED LOSITIVE INTEGER 
LRCM*i FOR FIRST CALL TO GETRAN 

SET LRCN GREATER THAN 1 FOR SUBSEQUENT CALLS TO GETRAN 

LACL-l FOR UNIFORM DISTRIBUTION 

LRCIXI2), DUM1 , DUM2 , NOT USED 1 

YFL “UNI FORML Y DISTRIBUTED RANDOM NUMBER 


COMPLEX ZAINR.l) 

COMMON /SE T/LRCL,LRCN, LRC 1X12) , 

IFIPCT ) 120,100,120 1 

100 IF(PCTB) 120,110,120 
110 IFIPHE) 120,140,120 
140 IF I 4MPL ) 120,145,120 

145 IFI AMPR ) 120,155,120 , , 

155 IFI AMP I 120.135,120 
120 P-PHE/57.296 
DO 130 J-1.N2 
K=(IC6*l)/2 
K1-IC6 

DO 130 I«l,Ni 

CALL GETRANILRCl Xil ) , LRCN.LRCL , DUM1 , YFL , 0JM2 ) 

IFILRCN.EQ. 1 )LRCN«LRCN*1 
E«2.0*P*( YFL-0.5) 

AA «ZA 1 1 , J ) 

AI»!0.,-1.)*ZA(I ,J| 

U« ABS I AR ) 

R2» ABSI A I ) 

A 3»R 1 »R2 

IF (R3.EQ.0.) GO TO 130 
R1-R1/R3 
R2-R2/R3 
A1«AA-E*AI 
A l *A I ♦E*AR 
AR«A1 

CALL GETRAN (LRC I X(1 ) , LRC N , LRCL , DUM1 , YFL , DJM2 ) 

E“i.0*2. J*PCT* (YFL-0.5) ♦PCTB 
AR»AR*c 
A I s A I *E 

CALL GETRANI LRC I XI 1 ) , LRC N , LRCL , DUM1 , YFL , DJM2 ) 

C AMP ERROR DIST IN PROPORTION TO RATIOS OF REAL AND IMAG TO I RE AL ♦ I MAG ) 
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IF (IC6.E0. 0 I GO TO 170 
I F I l-K ) 175,175»185 
175 AMP-AMPL 
GO TO 170 
185 AMP»AHPR 

1 73 t=2.*<tFL-.5 I »AMP*Ri 
AR«AR*£ 

CALL GETRANURCl X<1 ) , LRCN ,LRCL, DUHl, YFL , DUM2 > 
E»2.*(YFL-.5 l*ANP*R2 
AI-AI *E 

160 2A ( | « Jl > CMPLX (AR.A1 ) 

IF ( IC6.EQ. 0 ) GO TO 130 
IF. U.NE. Kl ) GO TO 130 
K*K»IC6 
K 1 a K 1 A I C6 
130 CONTINUE 
135 RETURN 
END 
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KllRC 3/13/73 2-DIM ORBITER 10,OOOLB/IN AT COORDINATES 5,13 
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* 
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DISPLACEMENT MOBILITY REAL* IMAGINARY FRE3« 92.00 HER TZ 
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KU.RC 3/1 3/73 2-01 1* ORB I TER 10.003LB/IN AT COORDINATES 5,13 

16 DEGREES OF FREEDOM 
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FORCE MATRIX AMP, PHASE IDEG) 
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CONSTRAINED RESPONSE MATRIX AMP, PHASE IOEG) 
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FREE BODY RESPONSE MATRIX AMP .PHASE (DEG) 
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